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§ ■ ABSTRACT 
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O . In recent years we have witnessed an explosion of photometric time-series 

data, collected for the purpose of finding a small number of rare sources, such as 
transiting extrasolar planets and gravitational microlenses. Once combed, these 
^ . data are often set aside, and are not further searched for the many other variable 

CL|. sources that they undoubtedly contain. To this end, we describe a pipeline that is 

Q . designed to systematically analyze such data, while requiring minimal user inter- 

^ . action. We ran our pipeline on a subset of the Trans-Atlantic Exoplanet Survey 

ci . dataset, and used it to identify and model 773 eclipsing binary systems. For each 

system we conducted a joint analysis of its light curve, colors, and theoretical 
, isochrones. This analysis provided us with estimates of the binary's absolute 

0^ . physical properties, including the masses and ages of their stellar components, 

00 ' as well as their physical separations and distances. We identified three types of 

O 
(N 



eclipsing binaries that are of particular interest and merit further observations. 
The first category includes 11 low-mass candidates, which may assist current ef- 
! forts to explain the discrepancies between the observation and the models of stars 

^ ' at the bottom of the main-sequence. The other two categories include 34 binaries 



with eccentric orbits, and 20 binaries with abnormal light curves. Finally, this 
uniform catalog enabled us to identify a number of relations that provide further 



5^ i constraints on binary population models and tidal circularization theory. 
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Introduction 



Since the mid 1990s there has been an explosion of large scale photometric variability 



survey s. The search for gravitational microlensing events, whi ch were predict e d by iPaczvnski 
fll986h. motiva t ed th e first wave of surveys [e.g. OGLE: ludalski et aP Jl994[) : EROS: 



Beauheu et al.l fll995h : DUO: lAlard fc GuibertI fll997h : MACHO: lAlcock et all fll998f )]. En- 
couraged by t h eir su ccess, additional surveys, searching for Gamma- Ray B ursts [e.g. ROTSE : 
Akerlof et al.l (120001 )] and general photometric variabilities [e.g. ASAS: iPojmanskil (119971 )] 
soon followed. 



S hortly thereafter, with the discovery of the first transiting extrasolar planet (jCharbonneau et al 



20001 : iHenrv et al 
e.g. OGLE-III: 



200G 



Mazeh et al.l boooh. a second wa ve of p hotometric surveys ensued 



^ _ _ Udalskil (l2003h : Tr ES:lAlonso et all (12004]): H AT: lBakos et al.l (120041): Super- 

WASP : [christian et 'allj2006h : XO: lMcCullough et al.l ^ 20061 ): for a review, see lCharbonneau et al 
( 120071 )]. Each of these projects involved intensive efforts to locate a few proverbial needles 



hidden in a very large data haystack. With few exceptions, once the "needles" were found, 
thus fulfilling the survey's original purpose, the many gigabytes of photometric light curves 
(LCs) collected were not made use of in any other way. In this paper we demonstrate how 
one can extract a great deal more information from these survey datasets, with comparably 
little additional effort, using automated pipelines. To this end, we have made all the software 
tools described in this paper freely available (see web links to the source code and working 
examples), and they are designed to be used with any LC dataset. 

In the upcomin g decade, a t hird w ave of ultra-la r ge gro und-based synoptic surveys 
[e.g. Pan-STARRS: IXaiser et al.l J2002h: L S ST: ITvsouI J2002h]. and ultra-sensitive space- 



based surveys [e.g. KEPLER : Borucki et al.l (Il997l ): COROT: iBadin fc The COROT Team 



(119981 ): GAIA: iGilmore et al.l (119981 )] are expected to come online. These surveys are designed 
to produce photometric datasets that will dwarf all preceding efforts. To make any efficient 
use of such large quantities of data, it will become imperative to have in place a large 
infrastructure of automated pipelines for performing even the most casual data mining query. 

In this paper, we focus exclusively on the identification and analysis of eclipsing bi- 
nary (EB) systems. EBs provide favorable targets, as they are abundant a nd ca n be well 
modeled using existing modeling programs [e.g. WD: IWilson &: Devinneyl (Il97ll ): EBOP: 
Popper &: Etzell (Il98ll )]. Once modeled, EBs can provide a wealth of useful astrophysical 
information, including constraints on binary component mass distributions, mass-radius- 
luminosity relations, and theories describing tidal circularization and synchronization. These 
findings, in turn, will likely have a direct impact on our understanding of star formation, stel- 
lar structure, and stellar dynamics. These physical distributions of close b inaries may even 
help solve open questions relating to the progenitors of Type la supernovae (llben fc Tutukov 
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1910 



1984 ) . In additional to these, EBs can be used as tools; both as distance indicators (IStebbing 



Paczvnskilll997f) and as sensitive detectors for tertiary companions via eclipse timing 



Deeg et aliboool : Iflolman fc Murravlboosl : kgol et al.l[2005h . 



In order to transform such large quantities of data into useful information, one must 
construct a robust and computationally efficient automated pipeline. Each step along the 
pipeline will either measure some property of the LC, or filter out LCs that do not belong, so 
as to reduce the congestion in the following, more computationally intensive steps. One can 
achieve substantial gains in speed by dividing the data into subsets, and processing them 
in parallel on multiple CPUs. The bottlenecks of the analysis are the steps that require 
user interaction. In our pipeline, we reduce user interaction to essentially yes/no decisions 
regarding the success of the EB models, and eliminate any need for interaction in all but two 
stages. We feel that this level of interaction provides good quality control, while minimizing 
its detrimental subjective effects. 

The data t hat we analy z ed or iginate from 10 fields of the Trans-atlantic Exoplanet 
Survey [TrES ; lAlonso et al.l (120041 )]. TrES employs a network of three automated tele- 
scopes to survey 6° x 6° fields-of-view. To avoid potential systematic noise we use the data 
from only one telescope , Sleuth, located at the Palomar Observatory in Southern California 
( 10 'Donovan et al.ll2004l ). This telescope has a 10 cm physical aperture and a photometric 
aperture of radius of 30". The number of LCs in each field ranges from 10405 to 26495 
(see Table [T]), for a total of 185445 LCs. The LCs consist of ~2000 r-band photometric 
measurements at a 9- minute cadence. These measurements were created by binning the 
image-subtraction results of 5 consecutive 90-second observations, thus improving their non- 
systematic photometric noise. As a result ~16% of the LCs have an RMS <1%, and ~38% of 
the LCs have an RMS <2% (see Table [2]). The calibration of TrES images, i dentification of 



stars therein, extraction, and decorrelation of the LC s is de scribed elsewhere ( jPunham et al. 



2004J : iMandushev et al.ll2005l : IO'Donovan et al.ll2006l . l2007l ) . TrES is currently an active sur- 
vey that is continuously observing new fields, though for this paper we have limited ourselves 
to these 10 fields. 



2. Method 



The pipeline we have developed is an extended version of the pipeline described by lDevor 

( 12005! ). At the heart of this analysis lie two computational routines that we have described 
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in earlier papers: the Detached Echpsing Binary Light curve fitteiQ [ DEBiL : IPevoil J2005h] 
and th e Method for Echpsing Component Identificatioij^ [MECI ; Devor &: Charbonneau 
fl2006al lb[)]. 3EBiL fits each LC to a geometric model of a detached EB (steps 3 and 5 below). 
This model consists of two luminous, limb-darkened spheres that orbit in a Newtonian 2- 
body orbit. MECI restricts the DEBiL fit along theoretical isochrones, and is thus able to 
create a physical model of each EB (step 9). This second model describes the masses and 
absolute magnitudes of the EB's stellar components, which are then used to determine the 
EB's distance and absolute separation. 



The pipeline consists of 10 steps. We elaborate on each of these steps below: 



1. Determine the period. 

2. If a distinct secondary eclipse is not observed, an entry with twice the period is added. 

3. Fit the orbital parameters with DEBiL. 

4. Fine-tune the period using eclipse timing. 

5. Refine the orbital parameters with DEBiL using the revised period. 

6. Remove contaminated LCs. 

7. Visually assess the quality of the EB models. 

8. Match the LC sources with external databases. 

9. Estimate the absolute physical properties of the binary components using MECI. 
10. Classify the resulting systems using both automatic and manual criteria. 



We use the same filtering criteria as described in iDevorl (120051 ). both for removing LCs 
that are not periodic (step 1) and then for removing non-EB LCs (step 3). Together, these 
automated filters remove approximately 97% of the input LCs. In addition to these filters, 
we perform stringent manual inspections (steps 7 and 10) whereby we removed all the LCs 
we were not confident were EBs. These inspections ultimately removed approximately 86% 



^The DEBiL source code, utilities, and running example files are available online at: 
ht t p : / / www ■ cf a. harvard . edu/ ~ j de vor/DEBiL . htiiil| 

^The MECI source code and running examples are available online at: 
http://www.cfa.harvard.edu/~jdevor/MECI.htmll 
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of the remaining LCs. Thus only 1 out of every 240 input LCs, were included in the final 
catalog. 



In step (1) we use both the Box-fitting Least Squares (BLS) peri od finder (IKovacs et al. 



1989 



2002 ) , and a version of the analysis of variances (AoV) period finder (jSchwarzenberg-Czerny 



19961 ) to identify the periodic LCs within the dataset and to measure their periods. 
In our AoV implementation, we scan periods from 0.1 days up to the duration of each LC. 
We then select the period that minimizes the varia nce of a linea r fit within 8 phase bins. 



We removed all systems with weak periodicities [see iDevorl (l2005l ) for details] , and with one 
exception (T-Lyr 1-14413), all the systems whose optimal period was found to be longer than 
half their LC duration. In this way we were able to filter out many of the non-periodic 
variables. 

The AoV algorithm is most effective in identifying the periods of LCs with long duration 
features, such as semi-detached EEs and pulsating stars. The BLS algorithm, in contrast, 
is effective at identifying periodic systems whose features span only a brief portion of the 
period, such as detached EEs and transiting planets (see Figure [1]). However, the BLS 
algorithm is easily fooled by outlier data points, identifying them as short duration features. 
For this reason the BLS algorithm has a significantly higher rate of false positives than AoV, 
especially for long periods, which have only a few cycles over the duration of the observations. 
Therefore we limit the search range of the BLS algorithm to periods shorter than 12 days, 
although as Figure [T] illustrates, its efficiency at locating EBs rapidly declines at periods 
greater than 10 days. 

In step (2), we address the ambiguity between EBs with identical components in a 
circular orbit, and EBs with extremely disparate components. The phased LC of EBs with 
identical components contains two identical eclipses, whereas the phased LC of EBs with 
disparate components will have a secondary eclipse below the photometric noise level. These 
two cases are degenerate, since doubling the period of a disparate system will result in a 
LC that looks like an equal-component system. In the pipeline, we handle this problem by 
doubling such entries; one with the period found in step (1), and another with twice that 
period. Both of these entries proceed through the pipeline independently. In many cases, 
after additional processing by the following steps, one of these entries will emerge as being far 
less likely than the other (see appendix [A]), at which point it is removed. But in cases where 
photometry alone cannot determine which is correct, one needs to perform spectroscopic 
follow up to break the ambiguity. In particular, a double-lined spectrum would support the 
equal-component hypothesis. 

Step (3) is performed using DEBiL, which fits the fractional radii (^1^2) and observed 
magnitudes {magi^2) of the EB's stellar components, their orbital inclination (z) and eccen- 
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tricity (e), and their epoch [to) and argument of periastron (u;). DEBiL first produces an 
initial guess for these parame ters, and then iteratively improves the fit using the downhill 
simp lex method (jNelder 1965 ) with simulated annealing ( Kirkpatrick et al. 1983 : Press et al. 



19921 ). 



In step (4), we fine-tune the period (P) using a method based on eclipse timingfl, which 
we describe below. In order to produce an accurate EB model in step (9) it is necessary to 
know the system's period with greater accuracy than that produced in step (1). If we neglect 
to fine-tune the period, the echpses may be out of phase with respect to one another, and so 
the phased eclipses will appear broadened. Our timing method employs the DEBiL model 
produced in step (3), and uses it to find the difference between the observed and calculated 
(O — C) eclipse epochs. This is done by minimizing the chi-squared fit of the model to the 
data points in each eclipse, while varying only the model's epoch of periastron. When the 
period estimate is off by a small quantity (AP), the O — C difference increases by AP each 
period. This change in the O — C over time can be measured from the slope of the linear 
regression, which is expected to equal AP/P. Thus measuring such an O — C slope will 
yield the desired period correction (see Figure [2]). 

If the EB has an eccentric orbit, the primary and secondary eclipse will separate on 
the O — C plot, and form two parallel lines with a vertical offset of At (see Figure [3]). 
We measure this offset and use it as a sensitive method to detect orbital eccentricities. In 
particular, the value of At constrains ecosu, which in turn provides a lower limit for the 
system's eccentricity (ITsesevichl Il973l ) : 



This formula assumes an orbital inclination of i = 90°, making it a good approximation 
for eclipsing binaries. We use this method, in combination with DEBiL, to identify the 
eccentric EBs in the catalog (see Table [3]). However, in cases where the eclipse timing 
measures |e costal < 0.005, or when the eccentricity is consistent with zero, we assume that 
the EB is non-eccentric, and model it using a circular orbit. We further discuss the physics 
of these systems in §3.2[ 

Step (5) is identical to step (3), except that it uses the revised period from step (4). 
This step provides an improved fit to the LCs, as evidenced by an improved chi-squared 
value in over 70% of the cases. 



■^The source code and running examples are available online at 
[http: / /www. cfa.harvard.edu/~jdevor/Timing. html 
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In step (6) we locate and remove non-EB sources that seem to be periodic due to 
photometric contamination by true EBs. Such contaminations result from overlapping point 
spread functions (PSF) that cause each source to partially blend into the other. These cases 
can be easily identified with a program that scans through pairs of target^ and selects the 
ones that both have similar periods (see description below) and are separated by an angle 
that is smaller than twice the PSF. We found 14 such pairs, all of which were separated by 
less than 41", which is well within twice the TrES PSF (60"), while the remaining pairs with 
similar periods were separated by over 450" . Upon inspection, all 14 of the pairs we found 
had similar eclipse shapes, indicating that we had no false positives. Among each pair, we 
identify the LC with shallower eclipses (in magnitudes) as being contaminated and remove 
it from the catalog. 

We define periods as being similar if the difference between them is smaller than their 
combined uncertainty. We estimate the period uncertainty using the relation: sp (x P'^ /T, 
where T is the time interval between the initial and the final observations. One arrives 
at this relation by noticing that when phasing the LC, the effect of any perturbation from 
the true period will grow linearly with the number of periods in the LC (see step 4). This 
amplified effect will become evident once it reaches some fraction of the period itself, in 
other words, when ep{T/P) oc P. A typical TrES LC with a revised period will have a 
proportionality constant of approximately 1/1000. In order to avoid missing contaminated 
pairs (false negatives) we adopt in this step, the extremely liberal proportionality constant 
of unity. 

In step (7) we conduct a visual inspection of all the LC fits. Most EBs were successfully 
modeled and were included into the catalog as is. About 1% of the LCs analyzed had 
misidentified periods, as a result of failures of the period finding method of step (1). In 
most of these cases the period finder indicated either a harmonic of the true period or a 
rational multiple of a solar or sidereal day. In such cases we use an interactive periodogra nfl 
to find the correct period and then reprocess the LCs through the pipeline. Some entries 
were misidentified at step (2) as being ambiguous, even though they have a detectable 
secondary eclipse or have slightly unequal eclipses. In these cases the erroneous doubled 
entry was removed. Lastly, some of the EBs were not fit sufficiently well with DEBiL in 
step (5). These cases were typically due to clustered outlier data points, systematic noise, or 
severe activity of a stellar component (e.g. flares or spots), which caused DEBiL to produce 



"^We ran a brute force scan, which required 0(-/V^) iterations. But by employing a data structure that 
can restrict the scan to nearby pairs, it is possible to perform this scan in only 0{N) iterations, assuming 
that such pairs are rare. 

■^LC, created by Grzegorz Pojmanski. 
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erroneous initial model parameters. These cases were typically handled by having DEBiL 
produce the initial model parameters from a more smoothed version of the LC. 

In step (8) we match each system, through its coo rdinates, with the corr esponding source 
in the Two Micron All Sky Survey catalog [2MASS ; (jSkrutskie et al.ll2006l )]. This was done 
to obtain both accurate target positions and observational magnitudes. These magnitude 
measurements are then used to derive the colors of each EB, which are incorporated into 
the MECI analysis, as well as to estimate the EB's distance modulus (step 9). To this 
end, 2MASS provides a unique combination of high astrometric accuracy (~0.1") together 
with high photometric accuracy (~0.015 mag) at multiple near-infrared bands, all while 
maintaining a decent photometric resolving power (~3"). By employing these near-infrared 
bands we both inherently reduce the detrimental effects of stellar reddening, and are able to 
correct for much of the remaining extinction by fitting for the Galactic interstellar absorbtion. 

In order to use the measurements from the 2MASS custom J, if, and Kg filters, we 
converted them to the equivalent ESO-filter values so that they could be compared to the 
isochrone table values used in the MECI anal ysis. This conversion was done using approx- 
imate linear transformations (jCarpenterl 120011 ) . However, the colors of three EBs (T-AndO- 
10336, T-Cygl-02304, and T-Perl-05205) were so anomalous that they did not permit a 
reasonable model solution, thus we chose not to include any color information in their MECI 
analyses. 

In addition to its brightness, we also look up each EB's proper motion. Although proper 
motion is not required for any of the pipeline analyses, it provides a useful verification for low- 
mass candidates (see §3.11) . These systems are expected to have large proper motions, since 
they must be nearby to be observable in this magnitude-limited survey. The most extreme 
such case in t he catalog is CM Dra conis (T-DraO-01363), which has a proper motion of over 
1300 mas/yr (ISalim fc Gouldll2003l ). and is probably the lowest mass system in our catalog. 
To this end, we match each system to the Second U.S. Naval Observatory CCD Astrograph 
Catalog [UCAC release 2.4 ; IZacharias et al.l (120041 )]. When there was no match with UCAC, 
we use the more comprehens ive but less accurat e U.S. Naval Observatory photographic sky 
survey [USNO-B release 1.0 : lMonet et al.l (120031 )]. These matches were made using the more 
accurate aforementioned adopted 2MASS coordinates. However, because of their increased 
observational depth, and the fact that some high-proper motion targets are expected to have 
moved multiple arcseconds in the intervening decades, we chose to match each target to the 
brightest (-R-band) source within 7.5". It should be noted that the position of CM Draconis 
shifted by more than 22" and had to be matched manually, though 90% of the matches were 
separated by less than 0.6", and 98% were separated by less than 2" (see Figure H]). 



The proper motions garnered from these databases can be combined with distance es- 
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timates [D), to calculate the absolute transverse velocity [vtr) of a given EB: 

vtr ^ 4.741 km/sec L ™ ) ( ) , (2) 
\ 1 mas /yr J \ 1 /epc / 

where PM is the system's angular proper motion. In the catalog we list the right 
ascension and declination components [PM^ and PMs, respectively), so as to allow one to 
compute the system's direction of motion in the sky. The value of PM can be computed 
from its components, using: PM^ = PM| + PM^ cos^ 6, where 6 is the system's declination. 
When applying this formula, one should be aware that USNO-B folds the cos 5 coefficient 
into its listed PMa, while UCAC does not. 

Finally, we incorporate the USNO-B photometric B- and P-magnitude measurements 
into our catalog to provide a rough estimate of the optical brightness of each target. USNO-B 
lists two independent measurements in each of these filter, however in some cases one or both 
of these measurements failed. When both measurements are available, we average them for 
improved accuracy. However, each measurement has a large photometric uncertainty of ~0.3 
mag, thus even these averaged values will have errors that are over an order of magnitude 
larger that the photometric measurements of 2MASS. For this reason, and because of the 
increased effect of stellar reddening, we chose not to incorporate these data into the MECI 
analysis. However, USNO-B's high photometric resolution (~r') enabled us to detect many 
sources that blended with our targets in the TrES exposures. By summing the i?-band fiuxes 
of all the USNO-B sources within 30" of each target, we estimated the fraction of third-light 
included in each LC (see Figure [5]). Note that this measure provides only a lower bound to 
the true third-light fraction, as some EBs are expected to have additional close hierarchical 
components that would not be resolved by USNO-B. For most of the catalog targets, the 
third-light fiux fraction was found to be small (<10%). We therefore conclude that stellar 
blending will usually have only a minor effect on the MECI analysis results, however users 
should be aware of the potential biases in the calculated properties of highly blended targets. 
Though it was not applied to this catalog, in principal, given a third-light flux fraction at a 
well-determined LC phase, one could correct for the effects of blending. 

In step (9) w e analyze the LCs with MECI. We refer the reader to the full description 
of this method in Devor fc Charbonneau J2006al lbl). and provide here only a brief outline. 



Given an observed EB LC and out-of-eclipse colors, MECI will iterate through a range of 
values for the EB age and the masses of its two components. By looking up their radii and 
luminosities in theoretical isochrone tables, MECI simulates the expected LC and combined 
colors, and selects the model that best matches the observations, as measured by the chi- 
squared statistic. Or more concisely, MECI searches the (Mi, M2, a(y'e)-parameter space 



- 10 - 



for the chi-squared global minimum of each EB. Figures [6] and [7] show constant-age slices 
through such a parameter space. Once found, the curvature of the global minimum along the 
parameter space axes is used to determine the uncertainties of the corresponding parameters. 

The MECI analysis makes two important assumptions. The first is tha t EB stellar com- 



pone nts are coeval, which has been shown to generally hold for close binaries (j Claret fc Willems 



2OO2I ). When this assumption is violated, MECI will often not be able to find an EB model 
that successfully reproduces the LC eclipses. Such systems, which may be of interest in 
their own right, make up ~3% of the catalog and are further discussed later in this section. 
The second assumption is that there is no significant reddening, or third-light blended into 
the observations (i.e. from a photometric binary or hierarchical triple). Such blending in 
the LC will make the eclipses shallower, which produces an effect very similar to that of 
the EB having a grazing orbit. Thus, it will cause the measured orbital inclination to be 
erroneous, although it should rarely otherwise affect the results of the MECI analysis signif- 
icantly. However, the MECI analysis is sensitive to color biases caused by stellar reddening 
and blending. 

We reduce both these biases by incorporating 2MASS colors (see step 8), which are both 
less suspectable to reddening than optical colors, and suffers from significantly less blending 
than TrES, as the radius of the 2MASS photometric aperture is ~20 times smaller that that 
of TrES. We then attempt to further mitigate this problem by an alyzing each EB twice, 



using d ifferent relative LC/color information weighting values [see iDevor &: Charbonneau 



( l2006bl ) for further details]. We first run MECI with the default weighting value {w = 10), 
and then run MECI again with an increased LC weighting {w = 100) thereby decreasing 
the relative color weighting. Finally, we adopt the solution that has a smaller reduced chi- 
square. Typically, the results of the two MECI analyses are very similar, indicating that the 
observed colors are consistent with the ones predicted by the theoretical isochrones. In such 
cases, the color information provides an important constraint, which significantly reduces 
the parameter uncertainties. However, when there is a significant color bias, the default 
model will not fit the observed data as well as the model that uses a reduced weighting of 
the color information. In such a case, the reduced color information model, which has a 
smaller chi-squared, is adopted. Following this procedure, we find that in ~9% of our EBs, 
the reduced color information model provided a better fit, indicating that while significant 
color-bias is uncommon, it is a source of error that should not be ignored. 



By default, we had MECI use the Yonsei-Yale flYi et al.ll200ll : iKim et al.ll2002f ) isochrone 
tables of solar metalhcity stars. Although they successfully describe stars in a wide range of 
masses, these tables become increasingly inaccurate for low-mass stars, as the stars become 
increasingly convective. For this reason we re-analyze EBs for which both components were 
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found to have masses below 0.75Mq, using instead the lBaraffe et al.l (119981 ) isochrone tables, 
assuming a convective mixing length equal to the pressure scale height. Our EB models also 
take into account the effects of the limb dark ening of each of the stellar compon ents. To 
this end we employ the ATLAS JKurucjll99^ and PHOENIX Jciaretlll998l . boooh tables of 
quadratic limb-darkening coefficients. 

As previously mentioned, o nce we know the absolute prop erties of an EB system, we 



are able to estimate its distance (IStebbing]ll910l: |Paczvnskilll997l ). and thus such systems can 



be considered standard candles. We use ICoxl (120001 ) for the extinction coefficients, assuming 
the standard Galactic ISM optical parameter, Ry = 3.1, to create the following system: 



magj-Magj = AM ag + 0.282 ■ A{V) (3) 
magH-Magu = AMag + 0.17Q ■ A{V) (4) 
magK-MagK = AMag + 0.108 ■ A{V) (5) 

Where AMag is the extinction-corrected distance modulus, and A{V) is the V-mag 
absorption due to Galactic interstellar extinction. The estimated distance can then be solved 
using: D = lOpc ■ lO'^^^"-9/^, Because we have three equations for only two unknowns, we 
adopt the solution that minimizes the sum of the squares of the residuals. In some cases 
we remove one of the bands as being an outlier (i.e. if it would have resulted in a negative 
absorption), after which we are still able to solve the systems. But in cases where we need 
to remove two bands, we set A{V) = in order to solve for the distance modulus. Although 
this method has a typical uncertainty of 10% to 20%, it can be applied to EBs that are far 
more distant and dim than are accessible in other methods, such as parallax measurement. 
It can be used to map broad features of the Galaxy, and identify binaries that are in the 
Galactic halo. This method can also be used on large group of systems, so that if the EBs are 
clustered, one can average their distances, and thus reduce the cluster distance uncertainty 
as the inverse square root of the number of systems. 

In step (10) we perform a final quality check for the EB model fits, and classify them 
into 7 groups: 

I. Eccentric: EBs with unequally-spaced eclipses 

II. Circular: EBs with equally-spaced but distinct eclipses 

III. Ambiguous-unequal: EBs with undetected secondary eclipses 

IV. Ambiguous-equal: EB with equally-spaced and indistinguishable echpses 
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V. Inverted: detached EBs that are not successfuUy modeled by MECI 
VI. Roche-lobe-filling: non-detached EBs that are filling at least one Roche-lobe 
VII. Abnormal: EBs with atypical out-of-eclipse distortions 

We list the model parameters for the EBs of groups I-IV in the electronic version of 
this catalog (see full description in appendix [B]). The EBs of groups V-VII could not be well 
modeled by MECI, therefore we list only their coordinates and periods, so that they can be 
followed- up. 

Figure [8] illustrates the period distribution of these seven groups. Note however that 
both the orbital geometry of EBs (eclipse probability oc P~^/^), and the limited duration of 
the TrES survey data (< 90 days ; varies from field to field ; see Tabled]), act to suppress the 
detection of binaries with longer periods. An added complication for single-telescope surveys 
is that about half of the EBs with periods close to an integer number of days will not be 
detectable, as they eclipse only durin g the daytini e. This EB distribution is consistent with 



the far deeper OGLE II field catalog (IDevorl 120051). where th e long tail of Roche-lobe-filling 



systems has recently been explained by iDerekas et al.l (120071 ) as being the result of a strong 



selection towards detecting eclipsing giant stars. 

Group [I] contains the eccentric EBs identified in step (4) as having centers of eclipse 
that are separated by a duration significantly different from half an orbital period (see 
Figures [HI [IHl and [TT]) . This criterion is sufficient for demonstrating eccentricity, but not 
necessary, since we miss systems for which cosu ~ (see equation [T]) . Fortunately, we are 
able to detect eccentricities in well-detached EBs with |ecosu;| > 0.005, using eclipse timing. 
Therefore, assuming that u is uniformly distributed, we are approximately 67% complete for 
e = 0.01, and over 92% complete for e = 0.04. In principle, it would be possible to be 100% 
complete for these systems by measuring the differences in their eclipse durations, however 



this measurement is known to be unreliable (lEtzell Il99ll ) and so would likely contaminate 



this group with false positives. Group [II] consists of all such circular-orbit EBs that were 
successfully fit by a single MECI model (see Figure [T2|) . 

EBs with only one detectable eclipse can potentially be modeled in two alternative ways. 
One way is to assume very unequal stellar components, which have a very shallow undetected 
secondary eclipse (group [III]). Since we cannot estimate the eccentricity of such systems, 
we assume that they have circular orbits. The other way is to assume that the period at 
hand is twice the correct value, and that the components are nearly equal (group [IV]). The 
entries of such ambiguous LCs were doubled in step (2), so that these two solutions would be 
independently processed through the pipeline (see Figure US]). Therefore, these two groups 
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have a one-to-one correspondence between them, although only one entry of each pair can be 
correct. Resolving this ambiguity may not always be possible without spectroscopic data. In 
some cases we were able to resolve this ambiguity using either a morphological or a physical 
approach. The morphological approach consists of manually examining the LCs of group [IV] 
for any asymmetries in the two eclipses (e.g. width, depth, or shape), or in the two plateaus 
between the eclipses (e.g. perturbations due to tidal effects, refections, or the "O'Connell 
effect"). The physical approach consists of applying our understanding of stellar evolution 
in order to exclude entries that cannot be explained through any coeval star pairing (see 
appendix Either way, once one of the two models has been eliminated, the other model 
is moved into group [II] and is adopted as a non-ambiguous solution. It is interesting to 
note that when analyzing the two models with MECI, the equal-component solution (group 
[IV]) has masses approximately equal to the primary component of the unequal-component 
solution (group [HI]). The mass of the unequal-component solution's secondary component 
will typically be the smallest value listed in the isochrone table, as this configuration will 
produce the least detectable secondary eclipse. 

Group [V] consists of detached EBs that cannot be modeled by two coeval stellar compo- 
nents. As mentioned earlier, we can reject the single-eclipse solution for EBs with sufficiently 
deep eclipses (see appendix This argument can be further extended to cases where we 
can detect both eclipses in the LC, but where one is far shallower than the other. In some 
cases, no two coeval main-sequence components will reproduce such a LC, but unlike the 
previous case, since both eclipses are seen, we cannot conclude that the period needs to be 
doubled. Such systems are likely to have had mass transfer from a sub-giant component onto 
a main-sequence component through Roche-lobe overflow, to the point where currently the 
main-sequ ence component has become significantly more massive and brighter than it was 



originally (jCrawfordI Il955l ) . This process will cause the components to effectively behave 
as non-coeval stars, even though they have in fact the same chronological age. In extreme 
cases, the originally lower-mass main-sequence component can become more massive than 
the sub-giant, and thus swap their original primary /secondary designations, so that the 
main-sequence component is now the primary component. We call such systems "inverted" 
EBs, and place them into group [V] (see Figure [T^ . This phenomenon is often referred to 
in the literature as the "Algol paradox", though we chose not to adopt this term so as to 
avoid confusing it with the term "Al gol- type EB" (EA), which i s defined by the General 



Catalogue of Variable Stars [GCVS ; (lKukarkinlll948l : ISamusll2006l )] as being the class of all 
well-detached EBs. 

Group [VI] contains the EBs that have at least one component filling its Roche-lobe 
(see Figure [T5l) . Such system cannot be well-fit by either DEBiL or MECI since they assume 
that the binary components are detached, and so neglect tidal and rotational distortions. 
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gravity darkening, and reflection effects. These systems must be separated from the rest of 
the catalog since their resulting best-fit models will be poor and therefore their evalu ated 
physical attributes will likely be erroneous. In a si milar fashi o n to iTamuz et al.l (120061 ) . we 



detect these systems automatically by applying the lEggletonl (Il983l ) approximation for the 



Roche-lobe radius, and place in group [VI] all the systems for which at least one of the EB 
components has filled its Roche-lobe (see Figure [16]), that is, if either one of the following 
two inequalities occurs: 



ri > 



r2 > 



0.49 g 



-2/3 



0.6g-2/3 + ln(l + g-i/3) 

0.49 g2/3 
0.6 g2/3 + In (1 + gi/3) ' 



or 



(6) 
(7) 



where q = M2/M1 is the EB components' mass ratio. Since we expect non-detached 
EBs to be biased towards evolving, higher mass stellar compo nents, we estimated q usin g 
the early-type mass-radius power law relation found in binaries (IGorda &: Svechnikovlll998l ): 
q ^ (r2/ri)^'^^^. Although in principle, we could have estimated q directly from the EB 
component masses resulting from the MECI analysis, we chose not to, since as stated above, 
the analysis of such systems is inaccurate. The analytic approximation we used, though 
crude, proved to be remarkably robust, as we found only 5 false negatives and no false 
positives when visually inspecting the LCs. We found ma ny r nore false po si tives/ negatives 
when using the alarm criteria suggested by iDevorl (120051 ) or iMazeh et al. J2006f ). both of 
which attempt to identify bad model fits by evaluating spatial correlations of the model's 
residuals. 

Finally, group [VII] contains systems visually identified as EBs (i.e. having LCs with 
periodic flux dips), yet having atypical LC perturbations that indicate the existence of ad- 
ditional physical phenomena (see Figures [T71 and [TH!) . For lack of a better descriptor, we 
call such systems "abnormal" (see further information in §3.3p . This group is different from 
the previous six in that we cannot automate their classification, and their selection is thus 
inherently subjective. In 15 of the 20 systems, we were able to approximately model the 
LCs, and included them in one of the aforementioned groups. In these cases users should 
be aware that these model may be biased by the phenomenon that brought about their LC 
distortion. 
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3. Results 

We identified and classified a total of 773 EB^. These systems consisted of 734 EBs with 
circular orbits, 34 detached EBs with eccentric orbits (group [1] ; Table [3]), and 5 unclassified 
abnormal EBs (group [Vll] ; Table [5]). We marked 15 of the detached EBs with circular orbits 
as also being abnormal. Of the 734 EBs with circular orbits, we classify 290 as unambiguous 
detached EBs (group [11] ; Table [7]), 103 as ambiguous detached EBs, for which we could 
not determine photometrically if they consisted of equal or disparate components (groups 
[III] and [IV] ; Table |6]), 23 as inverted EBs (group [V] ; Table [8]), and 318 as non-detached 
(group [VI] ; Table H]). With the exception of the abnormal EBs, which were selected by eye, 
we use an automated method to classify each of these groups (see §2] for details). Our mass 
estimates for the primary and secondary components are plotted in Figure [191 

The EB discovery yield (the fraction of LCs found to be EBs), varies greatly from field 
to field, ranging from 0.72% for Cygnus, to 0.15% for Corona Borealis (see Table [2]). This 
variation is strongly correlated with Galactic latitude, where fields near the Galactic plane 
have larger discovery yields than those that are farther from it (see Figure [20]) . This effect 
is likely due to the fact that fields closer to the Galactic plane contain a higher fraction of 
early-type stars. These early-type stars are both physically larger, making them more likely 
to be eclipsed, and are more luminous, which has them produce brighter and less noisy LCs, 
thereby enabling the detection of EBs with shallower eclipses. Furthermore, much of the 
residual scatter can be attributed to the variation in the observed duration of each field (see 
Table [1]). That is, we find additional EBs, with longer periods, in fields that were observed 
for a longer duration. 

Currently, 88 of the cataloged EBs (11%) appear in either the International Variable 
Star Index0 (VSX), or in the SIMBAd|^ astronomical database (Table [9]). However, only 49 
systems (6%) have been identified as being variable. Not surprisingly, with few exceptions, 
these targets were among the brightest sources of the catalog. Using only photometry, it 
is often notoriously difficult to distinguish non-detached EBs from pulsating variables that 
vary sinusoidally in time, such as type-C RR Lyrae. Furthermore, unevenly spotted stars 
may also cause false positive identifications, especially in surveys with shorter durations. 
Ultimately, spectroscopic follow-up will always be necessary to confirm the identification of 



^The observed LCs, fitted models, and model residuals of each of these EBs are shown at: 
ht t p : / / www . cf a. harvard . edu/ ~ j de vor / C at alog . h t ml] 

''Maintained by the American Association of Variable Star Observers (AAVSO). 

^Maintained by the Centre de Donnes astronomiques de Strasbom'g (CDS). 
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such variables. 

We highlight three groups of EBs as potentially having special importance as test beds 
for current theory. For more accurate properties, these EBs will likely need to be followed up 
both photometrically and spectroscopically. The brightness of these EBs will considerably 
facilitate their follow-up. 



3.1. Low- mass EBs 



The first group consists of 11 low-mass EB candidates, including 10 newly discovered 
EBs with either K or M-dwarf stellar components. Our criteria for selecting these bina- 
ries were that they be well-detached, and that both components have estimated masses 
below O.75M0 (see Table [TO] and Figure [2Tll . Currently, only 7 such detached l ow-mass EBs 



have been confirme d [YY Ger n (IKron 



1952 



Torres fc Ribas 



Metcalfe et al.lll996h. CU Cnc JPelfosse et alJ|l999l : lRibasH2003l ), T-HerO-07621 JCreevev et al 



20021 ). CM Dra 



Lacv 



1977b 



2005h . GU Boo jLopez-Morales fc Ribasl boosf). NSVS01031772 Jhopez-Morales et al.lbooeL 
and UNSW-TR-2 jYoung et al.ll2006h ]. 



Despit e a great deal of work that has been done to understand the structure of low-mass 
stars [e.g. Chabrier fc Baraff^ (200o[)], models continue to underest i mate their r adii by as 



much as 15% fcacy 



1977a 



Torres fc Ribasll2002l : ICreevev et al.ll2005l : lRibasll2006h . a signifi 



cant discrepancy conside ring that for solar-typ e stars the agreement with the observations is 
typically within 1 — 2% (jAndersenlll99ll . 119981 ). In recent years an intriguing hypothesis has 
been put forwar d that strong magnetic fields may have bloated these stars through chromo - 



spheric activi ty (RibasI 



Furthermore, Torres et al 



200' 



Torres et al.ll2006l : iLopez- Morales! 120071 : IChabrier et all 120071 1. 



J2006h ind that such bloating occurs even for stars with nearly a 



solar mass, and suggest that this effect may also be due to magnetically induced convective 
disruption. In either case, these radius discrepancies should diminish for widely separated 
binaries with long periods, as they become non-synchronous and thus rotate slower, which 
according to dynamo theory would reduce the strength of their magnetic fields. 

Unfortunately, the small number of well characterized low-mass EBs makes it difficult 
to provide strong observational constraints to theory. Despite the fact that such stars make 
up the majority of the Galactic stellar population, their intrinsic faintness renders them 
extremely rare objects in magnitude-limited surveys. In addition, once found, their low fiux 
severely limits the ability of observing their spectra with both sufficiently high resolution and 
a high signal-to-noise ratio. To this end, the fact that the TrES survey was made with small- 
aperture telescopes is a great advantage, as any low-mass EB candidate found is guaranteed 
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to be bright, and thus require only moderate-aperture telescopes for their follow-up. Thus 
we propose multi-epoch spectroscopic study of the systems listed here, in order to confirm 
their low mass and to estimate their physical properties with an accuracy sufficient to test 
models of stellar structure. Moreover, two of our candidates (T-Cygl-12664 and T-CasO- 
10450), if they are in fact ambiguous-equal (group [IV]), have periods greater than 8 days, 
making them prime targets for testing the aforementioned magnetic-bloating hypothesis. 



3.2. Eccentric EBs 



The second group of EBs consists of 34 binaries with eccentric orbits (see Table [31 and 
Figures [9| [T0| and fTTj) . We were able to reliably measure values of |ecosci;| as low as ~ 0.005 
by using the eclipse timing technique (see ^and Figure [3]). Since this measure provides a 
lower limit to the eccentricity, it is well-suited to identify eccentric EBs, even though the 
actual value of the eccentricity may be uncertain. As mentioned earlier, in an effort to 
avoid false-positives, we do not include in this group EBs whose eclipse timing measures 
lecoscijl < 0.005, or EBs with an eccentricity consistent with zero. 

Our interest in t hese eccentric binaries stems from their potential to constrain tidal 
circularization theory (jDarwinlll879l ). This theory describes how the eccentricity of a binary 
orbit decays over time due to tidal dissipation, with a characteristic timescale {tdrc) that is 
a function of the components' stellar structure and orbital separation. As long as the com- 
ponents' stellar structure remains unchanged, the orbital eccentricity is expected to decay 
approximately exponentially over time [e oc exp(— t/tcirc)]- Howeve r, once the component s 
evolve off the main sequence, this time scale may vary considerably (jZahn &: Bouchetlll989l ). 
Thus, to understand the circularization history of binaries with circularization timescales 
similar to or larger than their evolutionary timescales, one must integrate over the evolu- 
tionary tracks of both stellar components. 

T hree alterna te tidal dissipation mechanisms have be en proposed: dynamic al tides 

(IZahn||l975l . 119771 ). equihbrium tides (IZahnI 119771 : lHutlll98ll ). and hydrodynamics (ITassoul 
19881 ). Despite its long period of development , the inherent di f ficulty of observing tidal 
dissipation has prevented definitive conclusions. IZahn fc Bouchetl ( 119891 ) add a further com- 
plication by maintaining that most of the orbital circularization process takes place at the 
beginning of the Hayashi phase, and that the eccentricity of a binary should then remain 
nearly constant throughout its lifetime on the main sequence. 

Observational tests of these tidal circularization theories, whereby td rn is measured 



statistically in coeval stellar populations, have so far proved inconclusive. iNorth fc Zahn 
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( 120031 ) found that short-period binaries in both the Large and Small Magellanic Clouds 
seem to have been circularized in agreement with the theory of dynamical tides. How- 
ever, iMeibom &: Mathieul (120051 ) show that with the exception of the Hyades, the stars in 



the clusters that they observed were considerably more circularized than any of the known 
dissipation mechanisms would predict. Furthermore, they find with a high degree of cer- 
tainty, that older clusters are more circularized than younger ones, thereby contradicting the 
Hayashi phase circularization model. 

Encouraged by the statistical effect of circularization that can be seen in our c atalog 
(Figure [22]) . we further estimated tdrc for each of the eccentric systems as follows. IZahn 
(119771 . 119781 ) provides an estimate for the orbital circularization timescale due to turbulent 
dissipation in stars possessing a convective envelope, assuming that corotation has been 
achieved: 



tr 



21g(l + g)/c2 V L 



1/3 



where M, R, L are the star's mass, radius and luminosity, and /c2 is the apsidal motion 
constant of the star, which is determined by its internal structure and dynamics. 

More massive stars, which do not have a convective envelope but rathe r develop a 



radiative envelope, are thought to circularize their orbit using radiative damping (jZahnlll975 



Claret fc Cunhal 119971 ). This is a far slower mechanism, whose circularization timescale can 
be estimated by: 



tr 



21g(l + g)iV6^2 \GM 



R' 



1/2 



a \ 21/2 
R^ 



(9) 



where E2 is the tidal torque constant of the star, and G is the universal gravitational 
constant. We can greatly sir nplify these expressions by applying Kepler's law [a^ = GM{1 + 
q){P/27iY], and adopt the ICoxl (120001 ) power law approximations for the main-sequence 
mass-radius and mass- luminosity relations. For the convective envelope case, we adopt the 
late- type mass-radius relation (M < l.SM©), and for the radiative envelope case we adopt 
the early-type mass-radius relation (M > 1.3Mq), thus arriving at: 



~ 1 1370Myr {E2/l0-^)-\-\l + qf^^ (P/dayf (M/Mo)"''^^ , M > 1.3Mq ^ ' 
Determining the value of k2 and E2 is the most difficult part of this exercise, since their 
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values are a function of the detailed stru cture and dynamics of the given star, which in turn 
changes significantly as the star evolves Jciaret fc Cunhalll997l : Iciaret fc Willemsll2002h . In 
our c a lculat ion, w e estin a ate th e se values by interpolatin g published theoretical tables [k2: 
ZahnI (119941 ). E^. IZahnI (119751 ): I Claret &: Cunhal (Il997l )]. Since both stellar components 



contribute to the circularization process, the combined circularization timescale becomes 
tcirc = l/(^,^rc,i + i^ rirc.?^' whcre the subscr ipts 1 and 2 refer to the primary and secondary 
binary components (jClaret fc Cunhal 1 199 71 ). In Table [3l we list the combined circularization 
timescale for each of the eccentric EBs we identify. 

The value of tdrc for most of the eccentric systems (21 of 34) is larger than the Hubble 
time, indicating that no significant circularization is expected to have taken place since 
they settled on the main sequence. About a quarter of the eccentric systems (8 of 34) 
have a tdrc smaller than the Hubble time but larger than 1 Gyr. While circularization is 
underway, the fact that they are still eccentric is consistent with theoretical expectations. 
The remaining systems (5 of 34) all have tdrc < 1 Gyr, have periods less than 3.3 days, and 
unless they are extremely young, require an explanation for their eccentric orbits. Two of 
these EBs (T-TauO-02487 and T-TauO-03916) are located near the star forming regions of 
Taurus, supporting the hypothesis that they are indeed young. However, this hypothesis 
does not seem to be adequate for T-CasO-02603, which has a period of only 2.2 days and 
tdrc — 0.26 Gyr, while possessing a large eccentricity of e ~ 0.25. An alternative explanation 
is that some of these binaries were once further apart, having larger orbital periods, and thus 
larger circularization timescales. These systems may have been involved in a comparably 
recent interaction with a third star (a collision or near miss), or have been influenced by 
repeated resonant perturbations of a tertiary companion. 

Finally, we would like to draw the reader's attention to our shortest-period eccentric 
EB, T-CasO-00394, whose period is a mere 1.7 days. Notably, this system is entirely con- 
sistent with theory, since its mass falls in a precarious gap, where the stellar envelopes of 
its components are no longer convective, yet their radiative envelopes are not sufficiently 
extended to produce significant tidal drag (see Figure [23]) . 



3.3. Abnormal EBs 

The third group of EBs consists of 20 abnormal systems (see Table and Figures [T71 
and [18]). While possessing the distinctive characteristics of EBs, these LCs stood out during 
manual inspection for a variety of reasons. These systems underline the difficulty of fully 
automating any LC pipeline, as any such system will inevitably need to recognize atypical 
EBs that were not encountered before. 
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The LCs we listed can be loosely classified into groups according to the way they deviate 
from a simple EB model. A few cases exhibited pulsation- like fluctuations that were not syn- 
chronized with the EB period (shorter-period: T-DraO-00398, longer-period: T-Lyrl-00359, 
T- Per 1-00750). These fluctuations may be due either to the activity of an EB component, or 
to a third star whose light is blended with the binary. In principle, one can identify the active 
star by examining the amplitude of the fluctuations during the eclipses. If the fluctuations 
originate from one of the components, their observed amplitude will be reduced when the 
component is being eclipsed. In such a case, if the fluctuations are due to pulsations, they can 
further provide independent con straints to the stellar properties through astro-seismological 



models (IMkrtichian et al.ll2004l ). To identify such fluctuating EBs one must subt r act th e 



fitted EB model from the LC, and evaluate the residuals [e.g. iPilecki fc Szczygiell (120071 )]. 
When the fluctuation period is fixed, one can simply search the residual LC using a peri- 
odogram, as was done in step (1) of our pipeline (see §2]). However, when the fluctuation 
period varies (i.e. non- coherent), as in the aforementioned LCs, one must employ alternate 
methods, since simply phasing their LC will not produce any discernable structure. For LCs 
with long-period fluctuations one can directly search the residuals for time dependencies, 
while for LCs with short-period fluctuations one can search the residuals for non-Gaussian 
distributions. However, in practice these measurements will likely not be robust, as there 
are many instrumental effects that can produce false positives. Thus, we employ a search 
for auto-correlations in the residual time series, which overcomes most instrumental effects, 
while providing a reliable indicator for many types of pseudo-periodic fluctuations. 

The remaining systems had LC distortions that appear to be synchronized with the 
orbital period. The source of these fluctuations is likely due to long-lasting surface inho- 
mogeneities on one or both of the rotationally synchronized components. When the LC 
has brief periodic episodes of darkening (T-AndO-11476, T-CasO-13944, T-Cygl-07584, T- 
DraO-04520), they can usually be explained as stable star spots, but brief periodic episodes 
of brightening (T-AndO-04594, T-HerO-08091), which may indicate the presence of stable 
hot-spots, are more difficult to interpret. This phenomenon is especially puzzling in the 
aforementioned two cases, in which the brightening episodes are briefer than one would 
expect from a persistent surface feature and repeat at the middle of both plateaus. 

When the two plateaus of a LC are not flat, they are usually symmetric about the center 
of the eclipses. This is due to the physical mirror symmetry about the line intersecting the 
binary components' centers. When the axis of symmetry does not coincide with the cen- 
ter of eclipse (T-AndO-00920, T-Cygl-08866, T-DraO-03105, T-Lyr 1-07584, T-Lyrl-15595), 
a phenomenon we term "eclipse offset" , we conclude that this symmetry must somehow be 
broken. This may occur if the EB components are not rotationally synchronized, or have a 
substantial tidal lag. Another form of this asymmetry can appear as an amplitude difference 
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between the two LC plateaus (T-HerO-03497, T-Lyrl-13166, T-Per 1-08789, T-UMaO-03090). 
This phenomenon, which was originally called the "periastron effect" and has since been 
renamed the "O'C onnell effect", ha s been k n own for over a century, and has been exten- 
sively studied [e.g. lO'Connelll (Il95ll ): iMilond (Il986l )]. Classic h ypotheses su ggest an uneven 
distributi on of circumstellar material orbiting with the binary (jStruvdll948l ) or surrounding 
the stars jMergental ei]ll950h. either of which could induce a preferential H absorption on 
one side. iBinnendijk (JL96Qi) was the first of many to suggested that the this asymmetry is 
due to subluminous regions of the stellar surface (i.e. star spots). However, this explanation 
also requires the stars to be rotationally synchronized, and for the spots to be stable over 
the duration of the observations. Alternate models abound, including a hot spot on one side 
of a component brought about through mass transfer from the other component, persistent 
star spots created by an off-axis magnetic field, and c ircumstellar mate rial being captured by 
the components and heating one side of both stars (iLiu fc Yang||2003l ). As with many phe- 
nomena that have multiple possible models, the true answer may involve a combination of a 
number of these mechanisms, and will likely vary from system to system (jPavidge fc Milone 



19M). 



Finally, a few particularly unusual LCs (T-DraO-03105, T-Lyrl-05984) display a very 
large difference between their eclipse durations. Although a moderate difference could be 
explained by an eccentric orbit, such extreme eccentricities in systems with such short orbital 
periods (0.5 and 1.5 days) are highly unlikely. 



4. Conclusions 

We presented a catalog of 773 eclipsing binaries found in 10 fields of the TrES survey, 
identified and analyzed using an automated pipeline. We described the pipeline we used to 
identify and model them. The pipeline was designed to be mostly automated, with manual 
inspections taking place only once the vast majority of non-EB LCs had been automati- 
cally filtered out. At the final stage of the pipeline we classified the EEs into 7 groups: 
eccentric, circular, ambiguous-equal, ambiguous-unequal, inverted, Roche-lobe-filling, and 
abnormal. The former four groups were all successfully modeled with our model fitting pro- 
gram. However, the latter three groups possessed significant additional physical phenomena 
(tidal distortions, mass-transfer, and surface activity), which did not conform to the simple 
detached-EB model we employed. 

We highlighted three groups of binaries, which may be of particular interest and warrant 
follow-up observations. These groups are: low-mass EEs, EEs with eccentric orbits, and 
abnormal EEs. The low-mass EEs (both components < O.75M0) allow one to probe the 
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mass-radius relation at the bottom of the main-sequence. Only 7 such EBs have previously 
been confirmed, and the physical properties of many of them are inconsistent with current 
theoretical models. Our group of 10 new candidates will likely provide considerable additional 
constraints to the models, and the discovery of 2 long period systems could help confirm a 
recent hypothesis that this inconsistency is due to stellar magnetic activity. The eccentric- 
orbit EBs may help confirm and constrain tidal circularization theory, as many of them have 
comparably short circularization timescales. We demonstrated that, as one would predict 
from the theory, the shortest period systems fall within a narrow range of masses, in which 
their stellar envelopes cease to be convective yet their envelopes are not extended enough to 
produce significant tidal drag. The abnormal EBs seem to show a plethora of effects that are 
indicative of asymmetries, stellar activity, persistent hot and cold spots, and a host of other 
physical phenomena. Some of these systems may require dedicated study to be properly 
understood. 

In the future, as LC datasets continue to grow, it will become increasingly necessary 
to use such automated pipelines to identify rare and interesting targets. Such systematic 
searches promise a wealth of data that can be used to test and constrain theories in regions of 
their parameter-space that were previously inaccessible. Furthermore, even once the physics 
of "vanilla" systems has been solved, more complex cases will emerge to challenge us to 
achieve a better understanding of how stars form, evolve and interact. 
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A. Rejecting single-eclipse EB Models 

An EB LC comprised of a deep eclipse and a very shallow eclipse, can occur in one of two 
ways. Either the secondary component is luminous but extremely small (e.g. a white dwarf 
observed in UV), thus producing a shallow primary eclipse, or the secondary component is 
comparably large but extr emely dim, thus pro ducing a shallow secondary eclipse. The first 



case, though possible [e.g. iMaxted et al.l (120041 )]. is extremely rare, and will have a signature 



"flat bottom" to the eclipse. We have not encountered such a LC in our dataset. The second 
case will have a rounded eclipse bottom, due to the primary component's limb darkening. 
Assuming this latter contingency, in which the secondary component is dark in comparison 
to the primary component, we can place a lower bound to its radius (-R2): 



R2 >Ri^l- lO-o-^^^^ai , (Al) 

where Ri is the radius of the primary component, and Amagi is the magnitude depth 
of the primary eclipse. Thus, if the eclipse is very deep, the size of the secondary component 
must approach the size of the primary component. However, coeval short-period detached 
EBs with components of similar sizes yet desperate luminosities are expected to be very rare, 
assuming they follow normal stelar evolution. Therefore, if only one eclipse is detected, and 
it is both rounded and sufficiently deep, we may conclude that this configuration entry is 
likely to be incorrect, and that the correct configuration has double the orbital period and 
produces two equal eclipses. Only when we cannot apply such a period doubling solution 
(i.e. when the secondary eclipse is detectable), do we resort to questioning our assumption 
of normal stelar evolution (see classification group V, described in 



B. Description of the Catalog Fields 

Due to the large size of the catalog we were only able to list small excerpts of it in the 
body of this paper. Readers interested in viewing the catalog in its entirety can download 
it electronically. Note that although the catalog lists 773 unique systems, each of the 103 
ambiguous EBs appear in both possible configurations (see 0, raising the total number of 
catalog entries to 876. Below, we briefiy describe the catalog's 38 columns. The column 
units, if any, are listed in square brackets. 

1. Category - the EB's classification (see 0. 

2. Binary name - the EB's designation, which is composed of its TrES field (see Table 
[1]) and index. 
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3. a - the EB's right ascension (J2000). 

4. 6 - the EB's dechnation (J2000). 

5. Period [days] - the EB's orbital period. 

6. Period uncertainty [days] - the uncertainty in the EB's orbital period. 

7. Massi [Mq] - the mass of the EB's primary (more massive) component. 

8. Massi uncertainty [Mq] - the uncertainty in the primary component's mass. 

9. Mass2 [Mq] - the mass of the EB's secondary (less massive) component. 

10. Mass2 uncertainty [Mq] - the uncertainty in the secondary component's mass. 

11. Age [Gyr] - the age of the EB (assumed to be coeval). 

12. Age uncertainty [Gyr] - the uncertainty in the EB's age. 



13. Score - a weighted reduced of the MECI model fit [see iDevor fc Charbonneau 
(l2006bl ) for further details]. 



14. Isoch rone source - isochrone tables used [Y2: iKim et al.l (120021 ). or Baraffe: iBaraffe et al. 
Jl998h ]. 

15. Color weighting - the r elative weight {w) of the LC fit, compared to the color fit [see 



Devor fc Charbonneaul (l2006bl ) for further details]. 



16. PM source - t h e dat abase that provided the p roper motion measu rement [UCAC: 



Zacharias et all (120041 ) . USNO-B: Monet et al.l (120031 1. or Salim03: ISalim k Gould 



(|2003[ )]. 

17. PMa [mas/yr] - the right ascension component of the EB's proper motion. 

18. PMs [mas/yr] - the declination component of the EB's proper motion. 

19. Location error [arcsec] - the distance between our listed location (columns 3 and 4) 
and the location listed by the proper motion database. 

20. mags - the USNO-B 5-band observational magnitude of the EB (average of both 
magnitude measurements, if available). 

21. magu - the USNO-B i?-band observational magnitude of the EB (average of both 
magnitude measurements, if available). 
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22. Third — light fraction - the fraction of third-light flux (i?-band) blended into the LC 
(i.e. the flux within 30", excluding the target, divided by the total flux within 30"). 

23. magj - the 2MASS observational J-band magnitude of the EB, converted to ESO 
J-band. 

24. magn - the 2MASS observational i7-band magnitude of the EB, converted to ESO 
if-band. 

25. magx - the 2MASS observational i^'^-band magnitude of the EB, converted to ESO 
K-hand. 

26. Magj - the absolute ESO J-band magnitude of the EB listed in the isochrone tables. 

27. MagH - the absolute ESO i7-band magnitude of the EB listed in the isochrone tables. 

28. Magx - the absolute ESO i^"-band magnitude of the EB listed in the isochrone tables. 

29. Distance [pc] - the distance to the EB, as calculated from the extinction-corrected 
distance modulus. 

30. A(y) - the EB's V-mag absorption due to Galactic interstellar extinction (assuming 
Rv = 3.1). 

31. sm{i) - the sine of the EB's orbital inclination. 

32. |ecos(ct;)| - a robust lower limit for the EB's eccentricity (see equation [T]). 

33. Eccentricity - the orbital eccentricity of the EB. 

34. Eccentricity uncertainty - the uncertainty in the orbital eccentricity of the EB. 

35. Amagi - the r-band primary (deeper) eclipse depth in magnitudes. 

36. Epochi - the Heliocentric Julian date (HJD) at the center of a primary eclipse, minus 
2400000. 

37. Amag2 - the r-band secondary (shallower) eclipse depth in magnitudes. 

38. Epoch2 - the Heliocentric Julian date (HJD) at the center of a secondary echpse, minus 
2400000. 
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Note that the value of the uncertainties (columns 6, 8 10, 12, and 34), were calculated by 
measuring the curvature of the parameter-space contour, near its minimum. This method 
implicitly assumes a Gaussian distribution of the parameter likelihood. If the likelihood 
distribution not Gaussian, but rather has a flattened (boxy) distribution, then the computed 
uncertainty becomes large. In extreme cases the estimated formal uncertainty can be larger 
than the measurement itself. 
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Table 1. Observational parameters of the TrES fields 



Field 


Constellation 


a 

(J2000)* 


S 

(J2000) 


Galactic 
coordinates (l,b) 


Starting 
epoch (HJD) 


Ending 
epoch (HJD) 


Duration 
(days) 


AndO 


Andromeda 


01 09 30.1255 


+47 14 30.453 


(126.11, -015.52) 


2452878.9 


2452934.9 


56.0 


GasO 


Cassiopeia 


00 39 09.8941 


+49 21 16.519 


(120.88, -013.47) 


2453250.8 


2453304.6 


53.8 


CrBO 


Corona Borcalis 


16 01 02.6616 


+33 18 12.634 


(053.49, +048.92) 


2453493.8 


2453536.8 


43.0 


Cygl 


Cygnus 


20 01 21.5633 


+50 06 16.902 


(084.49, +010.28) 


2453170.7 


2453250.0 


79.3 


DraO 


Draco 


16 45 17.8177 


+56 46 54.686 


(085.68, +039.53) 


2453093.8 


2453163.0 


69.2 


HerO 


Hercules 


16 49 14.2186 


+45 58 69.963 


(071.61, +039.96) 


2452769.9 


2452822.0 


62.1 


Lyrl 


Lyra 


19 01 26.3713 


+46 56 05.325 


(077.15, +017.86) 


2453541.8 


2453616.7 


74.9 


Perl 


Perseus 


03 41 07.8581 


+37 34 48.712 


(156.37, -014.04) 


2453312.8 


2453402.8 


90.0 


TauO 


Taurus 


04 20 21.2157 


+27 21 02.713 


(169.83, -015.94) 


2453702.7 


2453770.9 


68.2 


UMaO 


Ursa Major 


09 52 06.3560 


+54 03 51.596 


(160.87, +047.70) 


2453402.9 


2453487.8 


84.9 



^ICRS 2000.0 coordinates of the guide star, which is located at the center of the field of view. 



Table 2. The number of sources and yield of the TrES fields 



Field 


Number 
of LCb 


Number of observations 
in each LC 


Fraction 
RMS < 1% 


Fraction 
RMS < 2% 


Found 
EBs 


EB discovery 
yield 


AndO 


26495 


2357 


16.5% 


40.4% 


111 


0.42%: 


CasO 


22615 


2069 


11.0% 


38.2% 


119 


0.53% 


CrBO 


18954 


1287 


11.0% 


22.4% 


28 


0.15% 


Cygl 


17439 


3256 


30.3% 


65.7% 


125 


0.72% 


DraO 


15227 


2000 


11.8% 


26.4% 


42 


0.28% 


HerO 


16916 


974 


16.8% 


35.0% 


28 


0.18% 


Lyrl 


22964 


2815 


19.4% 


49.0% 


135 


0.59% 


Perl 


20988 


1647 


15.9% 


38.4% 


93 


0.44% 


TauO 


14442 


1171 


13.1% 


32.5% 


68 


0.47% 


UMaO 


10405 


1343 


13.6% 


29.5% 


24 


0.23% 



Table 3. Eccentric EBs 



Object a (J2000) (5 (J2000) Period [daj/s]" | e cos ti) | ti,„i„g'' | e cos cj | ^^opt^^^ e"* Mi/Mq M^/Mq age [Gj/r] t^i^^[Gyr] 



T-AndO-04144 


01 


17 


35 


.247 


49 


46 


16. 


,97 


7.869 


0.0072 


0. 


,0068 


14+0"** 
^ — 0.08 


0.84 (-1) = 


0.54 (-1) 


10.0 (-3) 


140 


T-AndO-17158 


01 


10 


09 


.143 


48 


18 


19. 


,68 


11.415 


0.0182 


0. 


,0180 




1.03 (-1) 


0.92 (-1) 


10.0 (-3) 


370 


T-AndO-24609 


00 


58 


29, 


.826 


49 


25 


08. 


,88 


17.997 


0.0794 


0. 


,0799 


lQ+°-^° 
^■-L^_0.02 


1.22 ± 0.10 


1.10 ± 0.30 


5.4 ± 12.0 


6400 


T-CasO-00394 


00 


32 


51, 


.608 


49 


19 


39. 


,36 


1.746 


0.0235 


0. 


0242 


024+° "^ 
"■^ — 001 


1.46 ± 0.01 


1.44 ± 0.01 


3.4 ± 0.3 


260 


T-CasO-02603 


00 


47 


08, 


.610 


50 


37 


19. 


,32 


2.217 


0.2098 


0. 


2143 


25+"-^* 
^■^'^_0.04 


1.25 ± 0.01 


0.75 ± 0.04 


5.4 ± 5.0 


0.26 


T-CasO- 04534 


00 


31 


04 


.585 


51 


52 


10. 


,88 


6.909 


0.0057 


0. 


,0048 


014+° °^ 


1.17 ± 0.04 


0.96 ± 0.15 


6.4 ± 6.9 


29 


T-CasO- 04947 


00 


47 


10, 


.336 


50 


45 


12. 


,36 


3.285 


0.0845 


0. 


,0845 


10+"-°* 
^■-L^ — O 02 


1.04 (-1) 


0.86 (-1) 


10.0 (-3) 


0.53 


T-CasO-05165 


00 


43 


59, 


.256 


51 


14 


00. 


,07 


2.359 


0.0311 


0. 


0327 


15+"°** 


1.50 ± 0.21 


0.76 ± 0.17 


2.7 ± 2.9 


0.34 


T-CasO-07630 


00 


37 


23 


.347 


47 


19 


20. 


,68 


5.869 


0.0200 


0. 


,0298 


0.038_j,;ppg 


1.15 ± 0.12 


0.87 ± 0.34 


5.9 ± 9.7 


13 


T-Cygl-01364 


20 


09 


38 


.211 


49 


05 


08. 


,02 


12.233 


N/A 


0. 


,3254 


53+"-°* 


1.03 ± 0.18 


0.50 ± 0.09 


0.4 ± 1.2 


1100 


T-Cygl-01373 


19 


55 


44, 


.105 


52 


13 


34. 


,61 


4.436 


0.0059 


0. 


,0054 


010+°"^ 


0.97 (-1) 


0.82 (-1) 


10.0 (-3) 


3.0 


T-Cygl-01994 


20 


03 


03, 


.111 


52 


42 


04. 


17 


14.482 


N/A 


0. 


,0107 




1.80 (-1) 


1.06 (-1) 


0.20 (-2) 


2300 


T-Cygl-02304 


20 


02 


04 


.388 


47 


34 


14. 


,75 


5.596 


0.1549 


0. 


1529 


23+" '-° 
^■^^-0.08 


2.20 ± 1.28 


0.72 ± 0.41 


0.7 ± 4.8 


46 


T-Cygl-02624 


19 


59 


25 


.926 


52 


23 


59. 


,91 


11.608 


0.0172 


0. 


0172 


0.068t°°| 


2.11 ± 0.05 


1.52 ± 0.03 


0.3 ± 0.1 


10^ 


T-Cygl-06677 


20 


07 


25, 


.526 


52 


22 


00. 


,54 


6.512 


0.0077 


0. 


,0069 




1.54 ± 0.20 


1.31 ± 0.22 


1.6 ± 1.9 


10<* 


T-Cygl-07248 


19 


54 


45, 


.937 


50 


24 


05. 


,32 


6.058 


0.1674 


0. 


1681 


17+"'"' 
^■^ ' —0 001 


1.68 ± 0.01 


0.87 ± 0.20 


2.0 ± 2.1 


33 


T-Cygl-07297 


20 


10 


46, 


.910 


49 


09 


29. 


,42 


11.613 


0.3019 


0. 


,3010 


n ,0+0.08 

0'3»_o 08 


0.97 (-1) 


0.55 (-1) 


10.0 (-3) 


830 


T-Cygl-07584 


19 


58 


58 


.012 


47 


38 


19. 


,26 


4.925 


0.0074 


0. 


0074 


0.022+J5°i 


0.94 (-1) 


0.90 (-1) 


10.0 (-3) 


4.7 


T-Cygl-09934 


20 


10 


44, 


.209 


51 


07 


51. 


77 


4.549 


0.0505 


0. 


,0501 


11 + 006 
— 06 


1.35 ± 0.64 


0.94 ± 0.41 


3.5 ± 5.6 


5.6 


T-Cygl-15752 


20 


13 


52, 


.454 


50 


52 


23. 


12 


9.372 


0.2402 


0. 


2402 


35+0-OS 


1.31 ± 0.04 


1.05 ± 0.11 


3.6 ± 4.9 


230 


T-Lyrl-09931 


18 


59 


08, 


.441 


48 


36 


00. 


,04 


11.632 


0.2207 


0. 


2209 


9c: + 0.04 
^^•^^-0.03 


0.91 ± 0.09 


0.67 ± 0.08 


2.7 ± 3.3 


730 


T-Lyrl-13841 


19 


06 


26, 


.558 


48 


28 


47. 


,04 


6.640 


0.0362 


0. 


,0362 


075+°-" 
"■^'°-0.04 


1.01 ± 0.27 


1.01 ± 0.24 


8.7 ± 13.1 


19 


T-Lyrl-14413 


19 


03 


41 


.143 


47 


36 


55. 


,78 


39.861 


0.5922 


0. 


,6240 


"■°^-0.006 


1.08 ± 0.34 


0.96 ± 0.26 


6.4 ± 18.9 


10^ 


T-Lyrl-14508 


18 


57 


40 


.271 


48 


40 


51. 


,28 


8.050 


0.1861 


0. 


1862 




1.34 ± 0.28 


1.20 ± 0.78 


2.9 ± 8.4 


220 


T-Lyrl-22359 


19 


10 


54 


.290 


49 


26 


06. 


,95 


12.319 


0.1990 


0. 


1984 


„ ,0 + 0.05 
^■■^■^-0.05 


0.97 ± 0.48 


0.97 ± 0.46 


6.9 ± 29.3 


550 


T-Perl-00769 


03 


31 


43, 


.915 


36 


31 


52. 


,36 


3.648 


0.0248 


0. 


,0263 


0.055tS;g| 


1.06 ± 0.01 


0.65 ± 0.03 


7.6 ± 2.1 


1.4 


T-Perl-04218 


03 


35 


33, 


.667 


40 


00 


49. 


18 


4.070 


0.0072 


0. 


,0079 




0.94 (-1) 


0.72 (-1) 


10.0 (-3) 


2.4 


T-Perl-05205 


03 


34 


19 


.432 


39 


32 


44. 


,41 


8.472 


0.0558 


0. 


,0592 


095+° " 


2.22 ± 0.01 


1.17 ± 0.28 


0.9 ± 1.7 


210 


T-Perl-08252 


03 


52 


00 


.670 


40 


03 


47. 


,73 


4.457 


0.0656 


0. 


,0645 


065+0 '"' 


1.56 ± 0.01 


1.40 ± 0.34 


2.4 ± 2.5 


10^ 


T-Perl-11424 


03 


47 


56, 


.473 


37 


31 


31. 


,83 


4.247 


0.2403 


0. 


2404 


24+0 02 
"•^^-0.006 


1.01 (-1) 


0.82 (-1) 


10.0 (-3) 


2.3 


T-Perl-17327 


03 


40 


45, 


.644 


34 


47 


57. 


,26 


3.946 


0.0332 


0. 


,0305 


n 069+0 2S 


1.10 ± 0.02 


1.09 ± 0.09 


8.4 ± 16.4 


1.2 


T-TauO-02487 


04 


21 


55, 


.933 


25 


35 


49. 


,28 


2.826 


0.0125 


0. 


,0054 


014+0 005 
"■"-'*-0.005 


1.74 ± 0.07 


1.01 ± 0.08 


1.6 ± 0.7 


0.39 


T-TauO-03916 


04 


23 


37, 


.351 


25 


46 


36. 


,00 


3.217 


0.0713 


0. 


,0706 


071+0 02 
"■O'i-0.004 


1.18 ± 0.01 


1.15 ± 0.03 


6.0 ± 4.4 


0.56 


T-UMaO-01822 


09 


53 


37 


.710 


52 


45 


44. 


72 


9.551 


0.1502 


0. 


1503 


31 + 02 
"■3-^-0.02 


1.01 ± 0.04 


1.00 ± 0.05 


8.3 ± 4.8 


130 



^Thc full precision of the measured period is listed in the electronic version of the catalog, together with its uncertainty and the epoch of the center of eclipse (see appendix iBt . 

'^Measurements made using the eclipse timing of step (4). Although these values are approximations, they do not suffer from nearly as much numerical error as the DEBiL 
measurement, and are therefore usually accurate. "N/A" marks LCs for which there were too few eclipses to be able to apply the timing method. 

*^The adopted value is a combination of the values measured with the timing method and with DEBiL. 

"^The uncertainties of the eccentricities are non-Gaussian, since they have a strict lower bound (e > |ecosLt;|). We truncated the quoted lower uncertainties at this value, though 
even at this truncated value the real uncertainty is beyond la. 

^When the most likely model is at the edge of the parameter space, MECI is not able to bound the solution, and therefore cannot estimate the uncertainties. We mark (-3) when 
the upper limit was reached, (-2) when the lower limit was reached, and {-1) if one of the other parameters is at its limit. 



Table 4. EBs that fill at least one of their Roche-lobes (first 20) 



Object a (J2000) S (J2000) Period [days] 



T- 


■AndO- 


03774 


00 


59 01.029 


46 


47 


17. 


.08 


1. 


.362 


T- 


■AndO- 


■04813 


01 


16 37.880 


47 


33 


23 


.43 


0. 


.662 


T- 


■AndO- 


■05140 


01 


03 22.258 


44 


66 


24. 


.31 


0. 


.981 


T- 


■AndO- 


■05153 


01 


18 48.278 


49 


39 


36 


.86 


0. 


.492 


T- 


■AndO- 


■05343 


00 


52 55.122 


48 


01 


37. 


.68 


0. 


.824 


T- 


■AndO- 


07638 


01 


09 


27.871 


49 


20 


33, 


.81 


0, 


.403 


T- 


■AndO- 


■ 07892 


00 


56 


i,-). 


,-)f)7 


48 


39 


10 


73 





.380 


T- 


AndO- 


08330 


01 


19 


15. 


,949 


48 


00 


17, 


.45 


0, 


.630 


T- 


■AndO- 


■08652 


00 


66 


58 ■ 


,855 


49 


06 


06, 


.00 


0, 


.336 


T- 


■AndO- 


■09528 


01 


22 


09 ■ 


,328 


47 


14 


29, 


.86 


0, 


.918 


T- 


■AndO- 


■10071 


01 


14 


50 ■ 


,412 


49 


17 


46, 


.28 


0, 


.387 


T- 


■AndO- 


■10206 


00 


55 


55^ 


724 


49 


49 


46, 


.56 


0, 


.869 


T- 


■AndO- 


■10511 


01 


19 


16, 


,430 


47 


07 


46, 


.27 


0, 


.563 


T- 


■AndO- 


■10722 


01 


04 


03, 


8,")9 


48 


37 


13, 


.04 


1 


.062 


T- 


■AndO- 


11354 


01 


18 


05. 


168 


46 


10 


14, 


.66 


0, 


.331 


T- 


■AndO- 


■11476 


01 


07 


32. 


,106 


45 


66 


44, 


.93 


6, 


.380 


T- 


■AndO- 


■11599 


01 


09 


28. 


113 


46 


18 


24, 


.86 


0, 


.280 


T- 


■AndO- 


11617 


01 


07 


28, 


020 


4.5 


22 


40 


35 





503 


T- 


■AndO- 


12453 


01 


17 


12, 


aiG 


46 


42 


35 


.43 





.448 


T- 


■AndO- 


■12769 


00 


52 


58. 


164 


44 


44 


11, 


.26 


0, 


.325 



00 

I 



Table 5. Abnormal EBs 



Object 




(J2000) 


S {J2000) 


Period 
(days) 


Classified 
in catalog? 


In SIMBAD/VSXr 
(see table 1^ 


Notes 


T-AndO-00920 


01 


17 


30, 


,677 


47 


03 


31. 


.61 


24.073 






Large asymmetric reflection (0.1 mag) offset eclipse 


T-AndO-04594 


01 


16 


10, 


,713 


48 


52 


18. 


.97 


3.910 


yea 




TUted^ ratcaus (s ots?) 


T-AndO- 11476 


01 


07 


32 


, 106 


45 


55 


44 


.93 


6.380 


yes 






T-CasO- 13944 


00 


29 


48, 


,990 


50 


49 


54. 


.06 


1 .739 


yes 




IrUgular^ecUpse depths 


T-Cygl-07584 


19 


58 


58, 


,012 


47 


38 


19. 


.26 


4.925 


yes 


no 


Large persistent spot 


T-Cygl-08866 


20 


08 


36, 


,448 


49 


29 


35. 


.79 


2.876 


yes 


no 


Offset eclipse^ 


T-DraO-00398 


16 


57 


33, 


,875 


59 


31 


51. 


.98 


1.046 


yes 


yes 


Active (has 0.2 mag fluctuations with periods of a few hours) 


T-DraO-03105 


16 


23 


02, 


,558 


59 


27 


23. 


.44 


0.485 


no 


yes 


Unequal eclipses'^ / semi-detached 


T-DraO- 04520 


16 


49 


57, 


,960 


56 


26 


45. 


.56 


3.113 


yes 


no 


Tilted plateaus (spots?) 


T-HerO-03497 


16 


52 


28, 


,391 


44 


51 


29. 


.63 


7.853 


yes 


no 


Unequal plateaus'^ 


T-HerO-08091 


16 


51 


52, 


,608 


47 


01 


47. 


.98 


2.694 


yes 


no 


Offset eclipse 


T-Lyrl-00359 


19 


15 


33, 


,695 


44 


37 


01. 


.30 


1.062 


yes 


yes 


Large recurring spots (~ 0.05 mag) 


T-Lyrl-02800 


19 


08 


18, 


,809 


47 


12 


48. 


.16 


4.876 


no 


no 


Semi-detached / unequal plateaus (spots?) 


T-Lyrl-05984 


18 


53 


50, 


,481 


45 


33 


20. 


.90 


1.470 


no 


no 


Unequal eclipses^ / semi-detached 


T-Lyrl-08305 


18 


56 


43, 


.798 


48 


07 


02. 


.86 


14.081 


yes 


no 


Large asymmetric reflection (0.05 mag) ; offset eclipse 


T-Lyrl-13166 


19 


02 


28, 


.120 


46 


58 


57. 


.75 


0.310 


no 


no 


Unequal plateaus ; misshapen eclipse (persistent spot?) 


T-Lyrl-15595 


19 


06 


05, 


.267 


49 


04 


08. 


.95 


9.477 


yes 


no 


Offset eclipse 


T-Pcrl-00750 


03 


47 


45, 


.543 


35 


00 


37. 


.08 


1.929 


yes 


yes 


Spots / active 


T-Perl-08789 


03 


54 


33, 


.282 


39 


07 


41 


.53 


2.645 


yes 


no 


Tilted plateaus 


T-UMaO-03090 


10 


08 


52, 


.180 


52 


45 


52. 


.49 


0.538 


yes 


yes 


Unequal plateaus 



00 

^Even when the LC plateaus are not flat, due to tidal distortion or reflections, the system's mirror symmetry normally guarantees that the eclipses will occur during a plateau 
minimum or maximum. When, as in these cases, the eclipses are significantly offset from the plateau minima/maxima we can conclude that some mechanism, perhaps severe I 
tidal lag, is breaking the system's symmetry. 

Might not be an EB. This LC could be due to non-sinusoidal pulsations. 

*^Thc two LC plateaus between the eclipses, have a significantly different mean magnitude. This may be due to one or both components being tidally locked, and having a 
persistent spot or surface temperature variation at specific longitudes. 



Table 6. Ambiguous EBs (first 10) 



Vcr. 




Object 


a 


(J2000) 


S (,12000) 


Period, [days] 




M2 IMr;\ 


age [fryr] 


A 


T- 


■AndO-00657 


01 


06 


06, 


,159 


47 


31 


59, 


,37 


6.725 


2.50 (-1) 


0.74 (-1) 


0.20 (-2) 


B 


T- 


■AndO-006B7 


01 


06 


06. 


.169 


47 


31 


69. 


.37 


13.456 


1.92 (-1)"= 


1.92 (-1) 


0.20 (-2) 


A 


T- 


■AndO-01203 


01 


03 


34. 


.746 


48 


32 


39. 


.27 


3.505 


1.86 ± 0.09 


0.66 ± 0.10 


0.89 ± 0.83 


B 


T- 


■AndO-01203 


01 


03 


34. 


.746 


48 


32 


39. 


.27 


7.011 


1.90 ± 0.12 


0.66 ± 0.19 


0.80 ± 1.13 


A 


T- 


■AndO- 06017 


01 


12 


48, 


.217 


49 


68 


07, 


.16 


2.543 


1.40 ± 0.35 


0.62 ± 0.77 


3.49 it 4.28 


B 


T- 


-AndO-06017 


01 


12 


48 


,217 


49 


58 


07, 


,16 


5.085 


1.18 ± 0.71 


1.12 ± 0.85 


3.12 ± 11.15 


A 


T- 


■ AndO-06500 


01 


25 


56 


()S3 


49 


23 


31 


,74 


5.337 


0.97 ± 0.20 


0.49 ± 0.53 


7.71 ± 16.33 


B 


T- 


•AndO-06500 


01 


25 


56, 


,083 


49 


23 


31, 


,74 


10.674 


1.01 ± 0.30 


0.93 ± 0.45 


0.74 ± 1.60 


A 


T- 


■AndO-06680 


00 


65 


48, 


.163 


46 


02 


48, 


.67 


4.551 


1.16 ± 0.04 


0.61 ± 0.20 


6.09 ± 8.91 


B 


T- 


■AndO-06680 


00 


68 


48, 


.163 


45 


02 


48, 


.57 


9.104 


1.16 ± 0.09 


0.96 ± 0.29 


6.24 ± 10.78 


A 


T- 


-AndO-08053 


01 


13 


69. 


.402 


45 


61 


43. 


.43 


4.116 


1.14 (-1) 


0.40 (-2) 


6.00 (-1) 


B 


T- 


■AndO-08053 


01 


13 


69, 


.402 


46 


61 


43, 


.43 


8.231 


1.09 ± 0.55 


1.05 ± 0.64 


3.22 ± 16.37 


A 


T- 


-AndO-08417 


01 


01 


39, 


,041 


45 


03 


32, 


,98 


2.053 


1.01 (-1) 


0.47 (-1) 


10.00 (-3) 


B 


T- 


-And()-08417 


01 


01 


39, 


,041 


45 


03 


32, 


,98 


4.106 


1.01 (-1) 


0.90 (-1) 


10.00 (-3) 


A 


T- 


-AndO-09365 


01 


01 


00, 


,459 


45 


14 


24, 


,77 


1.887 


1.05 ± 0.03 


0.43 ± 0.39 


8.74 ± 16.06 


B 


T- 


■AndO-09365 


01 


01 


00, 


.469 


46 


14 


24, 


.77 


3.774 


1.05 ± 0.05 


0.93 ± 0.52 


9.47 ± 23.55 


A 


T- 


■AndO-10518 


01 


07 


44, 


,417 


48 


44 


68, 


.11 


0.194 


0.90 (-1) 


0.40 (-2) 


0.40 (-1) 


B 




■ AndO-10518 


01 


07 


44 


,417 


48 


44 


68, 


.11 


0.387 


0.45 ± 0.27 


0,45 ± 0,28 


0,27 ± 0.54 


A 




■ AndO- 11453 


01 


05 


42 


,744 


44 


54 


02, 


.26 


0.784 


1.12 (-1) 


0.40 (-2) 


7.00 (-1) 


B 


T- 


■AndO-11453 


01 


05 


42, 


.744 


44 


54 


02, 


.26 


1.568 


1.02 ± 0.43 


1.01 ± 0.32 


8.81 ± 14.54 



Unequal eclipse model, assuming an unseen secondary eclipse. 
^Equal eclipse model, with double the period of the unequal model. 

"^When the most likely model is at the edge of the parameter space, MECI is not able to bound the solution, and therefore 
cannot estimate the uncertainties. We mark (-3) when the upper limit was reached, (-2) when the lower limit was reached, 
and (-1) if one of the other parameter is at its limit. 



Table 7. Circular EBs (first 20) 



Object 


a 


(J2000) 


5 (J2000) 


Period [doj/s] 


Mi/Mq 


M2/MQ 


age [Gyr] 


Proper motion 
source catalog 


PMo, 
[MAS/year] 


[MAS/year] 


T-AndO- 00194 


01 


20 


12,816 


48 


36 


41.36 


2.145 


2.07 ± 0.02 


0.97 ± 0.02 


0.59 ± 0.12 


UCAC 


28.4 


-12.2 


T-AndO- 00459 


01 


11 


24.846 


46 


57 


49.44 


3.655 


1.20 ± 0.01 


1.19 ± 0.01 


5.35 ± 1.13 


UCAC 


-1.6 


-20.6 


T-AndO-00745 


01 


03 


45.076 


44 


50 


41.14 


2.851 


1.86 ± 0.23 


1.02 ± 0.22 


1.04 ± 0.72 


UCAC 


-6.6 


-4.8 


T-And()-()14fil 


01 


06 


15.353 


45 


08 


25.66 


5.613 


1.47 ± 0.01 


1.45 ± 0.08 


2.76 ± 2.72 


UCAC 


-11.4 


2.8 


T- And 0-01 554 


01 


17 


04.999 


45 


54 


06.20 


1.316 


0.90 (-1)'^ 


0.84 (-1) 


Xo 00 (-3) 


UCAC 


-44.6 


-40.8 


T-AndO-01597 


01 


10 


32.071 


46 


49 


53.18 


3.503 


1.55 ± 0.03 


1.54 ± 0.01 


2.37 ± 0.76 


i;cAC 


2.9 


-5.5 


T-AndO-02462 


01 


IS 


00.594 


49 


27 


12.47 


3.069 


1.97 ± 0.69 


1.10 ± 1.31 


1.02 ± 1.58 


UCAC 


5.8 


-1.1 


T-AndO-02699 


01 


06 


44.813 


47 


31 


08.61 


1.759 


1.18 ± 0.02 


0.63 ± 0.07 


5.21 ± 3.37 


UCAC 


0.2 


-6.8 


T-Anda-02798 


01 


21 


18.346 


48 


48 


05.63 


2.860 


1.04 ± 0.10 


0.66 ± 0.13 


6.14 ± 9.51 


UCAC 


6.3 


-8.1 


T-Anda-03526 


01 


20 


17.451 


47 


39 


23.32 


1.536 


1.04 ± 0.02 


0.84 ± 0.02 


6.29 ± 2.37 


UCAC 


17.9 


-11.1 


T-And0-O404fi 


00 


55 


20.157 


47 


44 


53.20 


3.916 


1.30 ± 0.09 


1.25 ± 0.12 


3.10 ± 4.31 


UCAC 


-3.8 


-7.3 


T-AndO-04694 


01 


16 


10.713 


48 


52 


18.97 


3.910 


1.05 (-1) 


0.82 (-1) 


10.00 (-3) 


UCAC 


1.5 


-1.9 


T-AndO- 04829 


01 


16 


16.228 


47 


45 


58.97 


0.678 


0.99 (-1) 


0.92 (-1) 


10.00 (-3) 


UCAC 


-23.8 


44.4 


T-Anda-05241 


00 


56 


34.679 


46 


37 


02.91 


1.454 


1.56 ± 0.01 


1.47 ± 0.31 


2.69 ± 7.01 


UCAC 


-4.6 


-0.5 


T-Anda-0537B 


01 


10 


68.226 


49 


52 


48.69 


1.640 


2.13 (-1) 


1.85 (-2) 


1.00 (-1) 


UCAC 


-6.3 


0.1 


T-Anda-05794 


01 


12 


11.763 


47 


32 


30.94 


1.053 


2.06 ± 0.19 


1.08 ± 0.64 


1.08 ± 1.92 


UCAC 


-0.4 


-1.4 


T-Anda-06039 


01 


23 


37.548 


48 


25 


37.73 


4.923 


1.22 ± 0.05 


1.08 ± 0.31 


5.33 ± 7.17 


UCAC 


-2.6 


-5.0 


T-Anda-06340 


01 


01 


55.269 


49 


18 


38.23 


5.437 


1.33 (-1) 


0.40 (-2) 


4.00 (-1) 


UCAC 


0.3 


-2.9 


T-Anda-06538 


01 


20 


58.907 


49 


29 


08.89 


18.669 


1.33 ± 0.15 


0.97 ± 0.17 


3.38 ± 3.45 


UCAC 


1.1 


-6.8 


T-AndO-06632 


01 


22 


36.840 


47 


52 


53.29 


1.669 


1.69 ± 0.01 


1.46 ± 0.24 


2.21 ± 1.02 


UCAC 


-7.2 


-7.6 



00 

^Whcn the most likely model is at the edge of the parameter space, MECI is not able to bound the solution, and therefore cannot estimate the uncertainties. We mark (-3) | 
when the upper limit was reached, (-2) when the lower limit was reached, and (-1) if one of the other parameter is at its limit. 



Table 8. Inverted EBs 



Object a (J2000) S (J2000) Period [days] 



T-AndO- 13663 


00 


69 


57, 


.881 


45 


03 


41. 


,63 


3. 


.342 


T-CasO-02069 


00 


49 


17 


.959 


50 


39 


02 


,92 


2 


,830 


T-CasO-03012 


00 


45 


41 


.832 


51 


01 


35 


,40 


1 


,108 


T-CasO- 04618 


00 


46 


22 


.661 


50 


39 


17, 


,57 


2 


,798 


T-Cas()-07780 


00 


34 


18 


.779 


52 


00 


35 


,72 


1 


,852 


T-CasO-19045 


00 


21 


44, 


.707 


50 


32 


29, 


,56 


0, 


,785 


T-CaBO-19668 


00 


48 


01, 


.342 


47 


06 


11, 


.68 


1, 


.848 


T-CasO-21651 


00 


26 


34, 


.895 


46 


38 


42, 


.69 


1, 


.166 


T-Cygl-01956 


19 


53 


29 


.106 


47 


48 


49 


,86 


2 


,045 


T-Cygl-02929 


20 


11 


57 


.009 


48 


07 


03 


,59 


4, 


,263 


T-CygH7342 


19 


49 


54 


.197 


50 


53 


28, 


,08 


2 


,220 


T-Hcr()-05469 


16 


54 


51 


.245 


43 


20 


35, 


,89 





,899 


T-Lyrl-04431 


19 


12 


16, 


.047 


49 


42 


23, 


,58 


0, 


,903 


T-Lyrl-06887 


18 


62 


10, 


.489 


47 


48 


16, 


.67 


1, 


.802 


T-Lyrl-07179 


18 


49 


14, 


.039 


45 


24 


38, 


.61 


1, 


.323 


T-Lyrl-10989 


19 


06 


22 


.791 


45 


41 


53, 


.82 


2 


,015 


T-Lyrl-11067 


18 


52 


53 


.489 


47 


51 


26, 


,58 


2 


,241 


T-Porl-04353 


03 


45 


04 


.887 


37 


47 


15 


91 


2 


953 


T-Pcrl-06993 


03 


40 


59 


.668 


39 


12 


35 


90 


2 


125 


T-Perl-09366 


03 


49 


20 


.305 


39 


55 


41 


,97 


2 


,374 


T-Perl- 12217 


03 


28 


59, 


.454 


37 


37 


42, 


.14 


1. 


.690 


T-TauO-00686 


04 


07 


13, 


.870 


29 


18 


32 


.44 


6, 


.361 


T-UMaO-00127 


09 


38 


06, 


.716 


56 


01 


07, 


.32 


0, 


.687 



CO 
CO 



Table 9. EBs that appear in either the VSX or the SIMBAD astronomical databases 



Category 


Object 




(J2000) 


5 (J2000) 


Spectral type 


Classification 


Identifiers 


Circular 


T-AndO-00194 


01 


20 


12, 


.816 


48 


36 


41, 


.36 


A5 


Star 


BD-t-47 378 ; GSC 03269-00662 ; SAO 37126 ; AG-|-48 143 
PPM 43886 ; TYC 3269-662-1 


Circular 


T-AndO-00459 


01 


11 


24, 


.845 


46 


57 


49, 


.44 


F8 


EB of Algol type 


CO And ; GSC 03268-00398 ; TYC 3268-398-1 ; BD-|-46 281 


Ambiguous 


T-AndO-00657 


01 


06 


06, 


.159 


47 


31 


59 


.37 


KO 


Star 


BD-h46 254 ; GSC 03267-01349 ; TYC 3267-1349-1 

AG+47 120 ; PPM 43637 


Circular 


T-And0-U0745 


01 


03 


43. 


,076 


44 


50 


41 


.14 




Star 


TYC 2811-470-1 ; GSC 02811-00470 


Ambiguous 


T-AndO-01203 


01 


03 


34, 


, 745 


48 


32 


39 


.27 




Star 


TYC 3267-1176-1 ; GSC 03267-01176 


Circular 


T-AndO-04046 


00 


55 


20, 


.157 


47 


44 


53, 


.20 




Star 


GPM 13.833991-1-47.748193 


Roche-fill 


T-AndO-05153 


01 


18 


48, 


.278 


49 


39 


36 


.86 




EB of W UMa type 


QW And 


Roche-fill 


T-AndO-05343 


00 


52 


55, 


.122 


48 


01 


37, 


.68 




Star 


GPM 13.232700-1-48.019757 


Ro chc-fill 


T-AndO-07S92 


00 


56 


15, 


.567 


48 


39 


10 


.73 




EB 


NSVS 3757820 


Circular 


T-AndO-23792 


00 


54 


09, 


.254 


47 


45 


19 


.91 




Star 


GPM 13.538629-1-47.755510 


Rochc-fill 


T-CasO-00170 


00 


53 


37, 


,847 


48 


43 


33 


.83 




Star 


TYC 3266-195-1 ; GSC 03266-00195 


Eccentric 


T-Cas()-()()394 


00 


32 


51, 


,608 


49 


19 


39 


.36 


B3 


EB of j3 Lyr type 


V381 Gas ; BD-j-48 lf)2 : BV 179 


Roche-fill 


T-Cas 0-00430 


00 


40 


06, 


.247 


50 


14 


15, 


.64 


K4 


EB of W UMa type 


V523 Cas ; GSC 03257-00167 ; WR 16 ; CSV 5867 
IRXS J004005. 0-1-501414 ; TYC 3257-167-1 


Circular 


T-CasO- 00640 


00 


47 


06, 


.277 


48 


31 


13 


.14 




Star 


TYC 3266-765-1 ; GSC 03266-00765 


Circular 


T-Cas 0-00792 


00 


48 


26, 


.554 


51 


35 


02. 


.52 




Star 


TYC 3274-664-1 ; GSC 03274-00664 


Roche-fill 


T-CasO-02013 


00 


4(J 


46, 


.427 


46 


56 


57 


.41 




Star 


TYC 3253-1767-1 ; GSC 03253-01767 


Inverted 


T-CasO-02069 


00 


49 


17, 


.959 


50 


39 


02, 


.92 




EB 


V385 Cas 


Roche-fill 


T-CasO-08802 


00 


51 


32, 


.351 


47 


16 


42, 


.57 




Star 


GPM 12.884787-1-47.278540 


Roche-fill 


T-CrBO-00654 


16 


00 


14, 


.507 


35 


12 


31 


.56 




EB of W UMa type 


AS CrB ; GSC 02579-01125 ; NSVS 7847829 
ROTSEl J160014.54-|-351228.4 


Roche-fill 


T-CrBO-00705 


15 


55 


51, 


.838 


33 


11 


00 


.39 




EB of W UMa type 


ROTSEl J155551. 87-1-331100. 5 


Roche-fill 


T-CrBO-01589 


16 


10 


09, 


.313 


35 


57 


30 


.57 




Variable of S Set type 


ROTSEl J161009. 33+355730.8 


Ro chc-fill 


T-CrBO-01605 


16 


00 


58, 


,472 


34 


18 


54 


.34 




EB of W UMa or RR Lyr-C 


NSVS 7848126 ; ROTSEl ,1160(358. 45-|-341854. 5 


Roche-fill 


T-CrBO-04254 


16 


09 


19, 


.589 


35 


32 


11, 


.48 




EB of W UMa type 


ROTSEl J160919. 62-1-353210. 8 


Circular 


T-Cygl-00246 


19 


44 


01, 


.777 


50 


13 


57, 


.42 




Star 


TYC 3565-643-1 ; GSC 03565-00643 


rtocne-iiii 


1 -Vjygl-UU^U^ 


19 


54 


39, 


.939 


50 


36 


41, 


.91 




Star 


1 OODD-oUO-l ; ^jtOLj UoiDOO-UUDUO 


Ambiguous 


T-Cygl-01385 


20 


15 


21, 


.936 


48 


17 


14, 


.14 




Star 


TYC 3576-2035-1 ; GSC 03576-02035 


Circular 


T-Cygl-01627 


19 


45 


20, 


.426 


51 


35 


07, 


.22 




star 


TYC 3569-1752-1 ; GSC 03569-01752 


Roche-fill 


T-Cygl-04652 


20 


07 


07, 


.305 


50 


34 


01, 


.34 




EB of W UMa type 


GSC 03567-01035 


Ro chc-fill 


T-Cygl-04852 


19 


51 


59, 


.208 


50 


05 


29 


.61 




EB of W UMa type 


NSVS 5645908 


Circular 


T-Cygl-09274 


20 


16 


06, 


.814 


51 


56 


26, 


.07 




EB of W UMa type 


V1189 Cyg ; CSV 8488 ; GSC 03584-01600 ; SON 7885 


Roche-fill 


T-Cygl-11279 


19 


59 


53, 


.377 


49 


23 


27, 


.86 




X-ray source 


IRXS J195954.0-K492318 


Roche-fill 


T-Cygl-12518 


19 


58 


15. 


.339 


48 


32 


15, 


.79 




Variable star 


Mis V1132 


Roche-fill 


T-Cygl-14514 


19 


48 


05. 


.077 


52 


51 


16, 


,25 




EB of W UMa or RR Lyr-C 


V997 Cyg ; GSC 03935-02233 ; ROTSEl J194804.79-F525117. 


Ambiguous 


T-DraO-00240 


17 


03 


52, 


.919 


57 


21 


55, 


.54 




Star 


TYC 3894-898-1 ; GSC 03894-00898 


Ambiguous 


T-DraO-00358 


16 


45 


38, 


.339 


54 


31 


32, 


.02 




Star 


TYC 3879-2689-1 ; GSC 03879-02689 


Circular 


T-DraO-00398 


16 


57 


33, 


.875 


59 


31 


51, 


.98 




EB of Algol type/X-ray source 


RX J1657. 5 + 5931 ; IRXS .7165733.5 + 593156 
VSX J165733.8+593151 ; GSC 03898-00272 


Roche-fill 


T-DraO-00405 


16 


27 


49, 


.103 


58 


50 


23, 


.30 




Star 


TYC 3884-1488-1 ; GSC 03884-01488 


Roche-fill 


T-DraO-00959 


16 


27 


44, 


.159 


56 


45 


59, 


.30 




EB of W UMa type/X-ray source 


NSVS 2827877 ; IRXS J162743.9+564557 


Circular 


T-DraO-01363 


16 


34 


20, 


.417 


57 


09 


48, 


.95 


M4.5V 


EB of BY Dra type 


CM Dra ; CSI+57-16335 1 ; LSPM J1634+5709 ; G 225-67 ; ' 



o 



T-DraO-01346 16 62 12.346 57 43 31.70 



High proper-motion Star 



EB of Algol type 



G 226-16 

IDS 16326+5721 A ; [RHG95] 2616 ; SBC7 580 ; CCDM J16343+5710A 
GJ 630.1 A ; LP 101-15 ; IDS 16325+5721 A ; [GKL99] 324 ; LHS 421 
2MASS .T16342040+5709439 ; CCABS 108 ; CABS 134 ; GEN# +9.80225067 
RX J1634.3+5709 ; IRXH J163421. 2+570941 ; IRXS J163421. 2+570933 
PM 16336+5715 ; USNO 168 ; USNO-Bl.O 1471-00307615 ; NLTT 43148 
BPS BS 16080-0095 ; VSX J165212. 3+674331 ; GSC 03886-00583 



Table 9 — Continued 



Category 


Object 


cc 


(J2000) 


.5 {J2000) 


Spectral type 


Classification 


Identifiers 


Roche-fill 


T-DraO-02224 


16 


30 


01 


.408 


54 


46 


66. 


.80 




Star 


BPS BS 16084-0159 


Circular 


T-DraO-03021 


17 


01 


03 


.618 


55 


14 


54, 


.70 




EB of Algol type 


VSX J170103. 5+551455 ; GSC 03890-01216 


Abnormal 


T-Dra()-t);-ilt)5 


16 


23 


02 


.558 


59 


27 


23 


.44 




X-ray source 


IRXS J162303. 6+592717 


Roche- fill 


T-DraO- 1)5259 


16 


41 


48 


.751 


56 


22 


34, 


.40 




EB of W UMa type 


VSX J164148. 7+562234 ; GSC 03882-02264 ; USNO-Bl.O 1463-0278621 


Ambiguous 


T-Hcr()-t)t)274 


17 


00 


51 


.150 


45 


25 


35 


.94 




Star 


TYC 3501-2245-1 ; GSC 03501-02245 


Roche-fill 


T-HerO-01086 


16 


48 


15 


.539 


44 


44 


28, 


.73 




EB of W UMa type 


GSC 03082-00896 ; NSVS 5252572 ; IRXS J164817. 3+444430 


Roche-fill 


T-HerO-03579 


16 


36 


47, 


.390 


45 


24 


68, 


.19 




EB of W UMa type 


GSC 03499-01631 


Inverted 


T-HerO-05469 


16 


64 


51 


.245 


43 


20 


36, 


.89 




EB 


V747 Her ; SVS 2066 


Circular 


T-Lyrl-00359 


19 


15 


33 


.695 


44 


37 


01, 


.30 


GOV 


EB 

X-ray source 


V2277 Cyg ; GSC 03133-01149 ; ROTSEl J191533. 92+443704.9 
BD+44 3087 ; ILFl+44 155 ; IRXS ,1191533.7+443704 


Circular 


T-Lyrl-()()687 


18 


55 


27 


.911 


47 


13 


41 


.76 




Star 


TYC 3544-1392-1 ; GSC 03544-01392 


Circular 


T-Lyrl-01013 


18 




03 


.963 


47 


49 


08 


39 




Star 


TYC 3544-2565-1 ; GSC 03544-02565 


Circular 


T-Lyrl-01439 


19 


06 


13, 


.439 


46 


57 


26 


.42 




Star 


TYC 3545-2716-1 ; GSC 03545-02716 


Circular 


T-Lyrl-02109 


18 


67 


35, 


.415 


45 


07 


44. 


.10 




Cepheid variable star 


ROTSEl J186736.99+450752.5 


Roche-fill 


T-Lyrl-02166 


19 


06 


07 


.448 


46 


16 


07. 


.61 




X-ray source 


IRXS J190604. 8+461512 


Roche-fill 


T-Lyrl-03173 


18 


59 


45 


.531 


47 


20 


07, 


.34 




EB of W UMa type 


ROTSEl ,1185945.43+472007.0 


Roche-fill 


T-Lyrl-03211 


IS 


45 


56 


.939 


47 


19 


09 


.54 




EB of W UMa type/X-ray source 


ROTSEl ,1184556,86+471914.4 ; IRXS J184567.9+471906 


Roche-fill 


T-Lyrl-()327() 


18 


57 


33 


.098 


48 


05 


22 


.49 




EB of W UMa type 


ROTSEl ,1185733,12+480522.5 


Roche- fill 


T-Lyr 1-03783 


18 


50 


12 


.684 


45 


35 


44, 


()5 






G P M 2 8 2 . 5 5 2 8 5 8 + 4 5 , 5 9 5 5 2 1 


Inverted 


T-Lyrl-04431 


19 


12 


16 


.047 


49 


42 


23 


.58 




EB of Algol type 


NSV 11822 ; GSC 03550-01770 ; NSVS 5678839 ; SON 9371 


Roche-fill 


T-Lyrl-06706 


18 


47 


57 


.211 


44 


38 


11, 


.30 




EB of W UMa type 


ROTSEl J184767. 18+443810, 8 


Inverted 


T-Lyrl-06887 


18 


62 


10 


.489 


47 


48 


16, 


.67 




EB of Algol type 


WX Dra ; AN 24.1925 


Rochc-flll 


T-Lyrl-06583 


18 


52 


26 


.837 


44 




20, 


.86 




EB 


ROTSEl J186226.63+445527,8 


Inverted 


T-Lyrl-(J7179 


18 


49 


14 


.t)39 


45 


24 


38 


61 




Star 


GPM 282,308484+45,410868 


Roche-fill 


T-Lyrl-08406 


18 


50 


06, 


.942 


45 


41 


05, 


.95 




Star 


GPM 282,528833+45,685035 


Roche-fill 


T-Lyr 1-10276 


18 


46 


55, 


.088 


45 


00 


52, 


.27 




EB of W UMa type 


V696 Lyr ; GPM 281.729421+46,014636 ; GSC 03640-00085 
ROTSEl J184664,98+450054,7 


Inverted 


T-Lyrl-10989 


19 


06 


22, 


.791 


45 


41 


53, 


.82 




EB of Algol type 


V612 Lyr ; SON 10931 


Roche-fill 


T-Lyrl- 11226 


18 


46 


21 


.748 


45 


63 


28, 


.79 




EB of W UMa type or « Set 


V694 Lyr ; GPM 281,340617+46,891326 ; GSC 03640-01842 
ROTSEl J184522, 47+455321,0 


Roche-fill 


T-Lyrl-12772 


18 


52 


25 


.096 


44 


55 


40, 


.23 




EB of W UMa type 


ROTSEl J185226, 53+445527,8 


Abnormal 


T-Lyrl-13166 


19 


02 


28, 


.120 


46 


58 


67, 


.76 


F9V 


EB 


V361 Lyr ; SON 9349 


Roche-fill 


T-Perl-00328 


03 


41 


57, 


.108 


39 


07 


29. 


.60 


G6 


EB of Algol type 


HD 276743 ; BD+38 787 ; GSC 02863-00766 ; TYC 2863-755-1 


Circular 


T-Per 1-00469 


03 


34 


57, 


.745 


39 


33 


18. 


.70 


G6 


Star 


HD 276647 ; GSC 02866-01996 ; TYC 2866-1996-1 


Circular 


T-Perl-00760 


03 


47 


46, 


.543 


35 


00 


37. 


.08 




Double or multiple star 


TYC 2364-2327-1 ; GSC 02364-02327 ; CCDM J03478+3601BC 
ADS 2771 BC ; BD+34 732B ; CSI+34 732 2 ; NSV 1302 


Roche-fill 


T-Pcrl-00974 


03 


34 


43 


.738 


38 


40 


22 


.22 


A 


Star 


HD 275481 


Circular 


T-Pcrl-()1218 


03 


42 


33 


.165 


39 


06 


03, 


.63 


A 


EB 


HU Per ; HD 276742 ; SVS 922 


Roche-fill 


T-Perl-01482 


03 


48 


45, 


.999 


35 


14 


10, 


.06 


FO 


Star 


HD 279025 


Circular 


T-Per 1-02697 


03 


44 


32, 


.202 


39 


69 


34. 


.94 


K4V 


T Tau type Star 


[LH98] 94 i IRXS J034432, 1+396937 ; ISWASP J034433, 95+396948,0 


Inverted 


T-Perl-04363 


03 


46 


04, 


.887 


37 


47 


16, 


.91 




EB of Algol type 


HV Per ; SVS 368 ; P 107 


Roche-fill 


T-TauO-00397 


04 


30 


09, 


.466 


25 


32 


27. 


.06 


A3 


EB of 9 Lyr type 


GW Tau ; SVS 1421 ; HD 283709 ; ASAS 043009+2632.4 


Inverted 


T-TauO-00686 


04 


07 


13, 


.870 


29 


18 


32, 


.44 




EB of Algol type 


IL Tau ; SON 9543 


Roche-fill 


T-TauO-00781 


04 


12 


51 


.218 


24 


41 


44, 


.26 


G9 


Eruptivc/T Tau-typc Star 


V1198 Tau ; NPM2+24,0013 ; IRXS J041250. 9+244201 
GSC 01819-00498 ; RX ,10412.8+2442 ; [WKS96] 14 


Roche-fin 


T-TauO-01262 


04 


16 


28, 


.109 


28 


07 


35, 


.81 


K7V 


Variable Star of Orion Type 


V1068 Tau ; EM StHA 25 ; JH 165 ; EM LkCa 4 
HBC 370 ; ASAS 041628+2807,6 


Roche-fill 


T-Tau0-01715 


04 


19 


26, 


.260 


28 


26 


14, 


.30 


K7V 


T Tau-type Star/X-ray source 


V819 Tau ; HBC 378 ; NAME WK X-Ray 1 ; IE 0416,3+2830 
IRAS C04162+2819 ; TAP 27 ; [MWF83] PI ; WK81 1 
IRXS J041926,l+282612 ; X 04163+283 


Roche-fill 


T-TauO-06463 


04 


07 


27, 


.415 


27 


61 


06, 


.36 




EB of W UMa type 


V1022 Tau ; HV 6199 ; NSV 1464 



Table 9 — Continued 



Category 


Object 


a (J2000) 


6 {J2000) 


Spectral type 


Classification 


Identifiers 


Inverted 


T-UMaO-00127 


09 38 06.716 


56 01 07.32 


A2V 


EB of Algol type 


VV UMa ; GEN# -1-0.05601395 ; HIP 47279 ; TYC 3810-1290-1 














CSC 03810-01290 ; SBC7 384 ; GCRV 6211 ; BD-l-56 1395 














HIC 47279 ; SVS 770 ; AAVSO 0931+56 


Circular 


T-UMaO-()()222 


10 07 18.023 


56 12 37.12 


AO 


Star 


HD 237866 ; GSC 03818-00504 ; SAO 27524 ; AG+56 778 ; HIC 49581 














BD+56 1432 ; HIP 49581 ; YZ 56 6209 ; TYC 3818-504-1 


Roche-fill 


T-UMaO-01701 


10 03 02.856 


55 47 53.34 




X-ray source 


RX J100303. 4-1-554752 ; [PTVOS] H22 ; [PTV98] P29 


Circular 


T-UMaO-03090 


10 08 52.180 


52 45 62.49 


K2e 


Star 


GSC 03815-01161 ; RIXOS 229-302 ; RX J100851.6-H624663 


Roche-fill 


T-UMaO-03108 


10 04 16.780 


54 12 02.83 




BB of W UMa type 


NSVS 2632137 



Table 10. Low-mass EB candidates (Mi^2 < O.75M0 ; sorted by mass) 



Category 


Object 


a 


(J2000) 


S {,12000) 


Period [days] 


Mi/Mq 




ago [Gyr] 


Proper motion 
source catalog^ 


mas/year 


PMs 
mas/year 


Circular 


T-DraO-01363'' 


16 


34 20.417 


57 09 48.95 


1.268 


0.27 ± 0.02 


0.24 ± 0.03 


1.6 ± 1.8 


Salim & Gould (20031 


-1121 


1186 


AmbigEq'^ 


T-AndO-10518 


01 


07 44.417 


48 44 58.11 


0.387 


0.45 ± 0.27 


0.45 ± 0.28 


0.3 ± 0.5 


UCAC 


2.7 


-2.0 


AmbigEq 


T-Cygl-12664 


19 


51 39.824 


48 19 55.38 


8.257 


0.50 ± 0.20 


0.48 ± 0.19 


0.3 ± 0.4 


USNO-B 


-18 


-8 


AmbigEq 


T-CrBO-14232 


16 


10 22.495 


33 57 52.33 


0.971 


0.60 ± 0.24 


0.55 ± 0.29 


4.4 ± 8.8 


UCAC 


-15.2 


-24.2 


AmbigEq 


T-CrBO-14543 


15 


57 45.926 


33 56 07.28 


1.506 


0.60 (-l)"* 


0.60 (-1) 


0.2 (-2) 


UCAC 


-13.9 


13.3 


Circular 


T-Pcrl-13685 


03 


53 51.217 


37 03 16.73 


0.384 


0.60 (-1) 


0.50 (-1) 


10.0 (-3) 


UCAC 


-24.1 


-15.9 


AmbigEq 


T-CrBO-10759 


15 


52 18.455 


30 35 32.13 


1.901 


0.63 ± 0.24 


0.82 ± 0.21 


7.3 ± 49.8 


UCAC 


3.6 


-19.4 


AmbigEq 
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^Where possible, wc used the more accurate UCAC catalog, otherwise wc reverted to the USNO-B catalog. Since they arc dim and nearby, wc expect most of the low-mass binaries to 
have comparably large proper motions. 

^T his binary is CM Draconis, which has been extensively studied and found to have a masses of Mi ^ 0.2307 ± O.OOIOM© and M2 ^ 0.2136 ± O.OOIOMq iLacvlflOTTbl : iMetcalfe et al l 
For consistency, we listed the MECI results, which are off by less than O.O4Af0 1.5cr). Wc also adopted an alternate proper motion estimate, as its USNO-B values seems to be 
erroneous, probably due to its very high angular velocity. 

'^For clarity we list for the ambiguous systems, only the solution with approximately equal components. But it is likely that at least a few of the ambiguous systems may be unequal, 
with half the period. Such cases can be identified as single-line spectroscopic binaries, with the secondary component being no larger than a few O.IA/q. 

'^When the most likely model is at the edge of the parameter space, MECI is not able to bound the solution, and therefore cannot estimate the uncertainties. We mark (-3) when the 
upper limit was reached, (-2) when the lower limit was reached, and (-1) if one of the other parameter is at its limit. 
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Fig. 1. — The fraction of the EBs in the final catalog found using the BLS algorithm and 
the AoV algorithm. The number of EBs in each bin is shown in Figure [HI The BLS method 
excelled at identifying EBs with short-duration eclipses (compared to the orbital period), 
which predominately occur at periods > 0.75 days. The AoV method fared better with EBs 
that have long-duration eclipses, which predominately occur in sub-day periods. The AoV 
method also does we ll with EBs with peri od longer than 10 days, which may be dominated by 
giant-giant binaries (jPerekas et al.ll2007l ). and so also have broad eclipses. This plot demon- 
strates the importance of using multiple independent methods of identifying EB, otherwise 
the results will have a significant selection effect that may bias any statistical results. 
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Fig. 2. — An eclipse timing plot produced in step (4), showing the O — C residuals of 
the primary eclipses (circles) and the secondary eclipses (triangles). Here, T-Lyrl-14962 is 
shown with an assumed period of 5.710660 days, as measured with an AoV periodogram. 
The slope of the residuals indicates that the assumed period is inaccurate. The grey solid 
line is predicted by the best circular-orbit model, whereas the dashed lines are predicted by 
the best eccentric-orbit model (compare to Figure [3]). After correction, we get a fine-tuned 
period of 5.712516(ia?/s. This 0.03% correction is small but significant in that without having 
had this correction, the eclipses would have smeared out and widened. 
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Fig. 3. — An eclipse timing plot for T-Cyg 1-01373, with an assumed period of 4.436013 (iays. 
In contrast to Figure [2], the slope here is consistent with zero, thus indicating that the period 
does not need to be fine-tuned. However, the O — C offset between the primary (circles) 
and secondary (triangle) eclipses is significant (1449 seconds), indicating that this EB has 
an eccentric orbit. The reduced chi-squared of the best circular-orbit model (grey solid line) 
is = 12.9, while the reduced chi-squared of the best eccentric-orbit model (dashed lines) 
is = 0.95. Applying the O — C timing offset to equation [1] provides a lower limit to the 
binary's orbital eccentricity: e > |ecosc(j| ~ 0.00594. 
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Fig. 4. — The distribution of the catalog position errors when matching targets to the proper 
motion databases. In some cases, the position errors are dominated by the motion of the 
EB during the intervening years. 
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Fig. 5. — The distribution of the i?-band third-hght flux fraction in the catalog LCs. This 
fraction was calculated by summing the fluxes of all the USNO-B sources within 30" of the 
target, excluding the target, and dividing this value by the total flux within 30", including 
the target. The resulting fraction ranges from to 1. 
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T-AndO-00745 (oge=1.0 Gyr; CIrculor) T-UMoO-01822 (age = 8,0 Gyr; Eccentric) 




Primary mass [Mi/M©] Primary mass [M1/M0] 

Fig. 6. — MECI likelihood contour plots of a typical circular-orbit EB (T-AndO-00745 ; left) 
and eccentric-orbit EB (T-UMaO-01822 ; right). There is no significant difference in the 
way MECI handles these cases, and both usually have a single contour minimum. The plots 
shown here have the ages set to the values that produced the lowest MECI minima. 




Fig. 7.— MECI likehhood contour plots of a typical ambiguous EB (T-TauO-03579). These 
plots show the effect of assuming that the binary components are equal (left) or unequal 
(right). Note that the equal-component solution will have a nearly symmetric contour around 
the diagonal, while the unequal-component solution can provide only an upper limit to the 
secondary component's mass, in this case M2 < 1 Mq. The plots shown here have the ages 
set to the values that produced the lowest MECI minima. 
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Fig. 8. — The EB orbital period distribution within the catalog. Each bin is subdivided 
to show the number of binaries belonging to each of the classification groups described in 
^ Note that the ambiguous-equal and ambiguous-unequal entries represent the same stars, 
with entries in the former group having double the period of the latter. Notice also how the 
Roche-lobe-filling EBs dominate the sub-day bins, and have a long tail stretching well above 
10-day periods. Furthermore, the circular orbit EBs have a period distribution peak of at 
~2 days, while the eccentric orbit EBs peak at ~5 days. This is likely due to the orbital 
circularization that occurs preferentially in short-period systems (see also Figures [22] and 
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Fig. 9. — Eccentric EBs (panel 1). Note how the secondary echpse is not at phase 0.5, as 
would be in circular orbit EBs. 
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Fig. 10. — Eccentric EBs (panel 2). 
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Fig. 11. — Eccentric EBs (panel 3). 
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Fig. 12. — Examples of unambiguous EBs with circular orbits, with their best-fit MECI 
models (solid line). 
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Fig. 13. — Examples of ambiguous EBs. Left column: assuming very unequal components. 
Right column: assuming approximately equal components with double the period. 
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Fig. 14. — Examples of EBs classified as inverted EBs. We included the unsuccessful best-fit 
MECI model (solid curve) as an approximate reference to illustrate the LC of a corresponding 
binary that has had no mass transfer. Notice how the model LC is unable to achieve a 
sufficiently deep primary eclipse, while producing a secondary eclipse that is too deep. 
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Fig. 15. — Examples of EBs that are assumed to have filled at least one of their Roche 
lobes. We included, for illustration purposes only, their best-fit MECI models (solid fine). 
These models were not adopted since they neglect tidal distortions, refiections, and gravity 
darkening effects, and so produce a poor fit to the data. 
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Fig. 16. — The criterion applied in equations [6] and [7] to determine whether one or both the 
EB components have filled their Roche lobe, and thus need to be placed into group (6). 
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Fig. 17. — LCs of abnormal EBs (panel 1). 
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Fig. 18. — LCs of abnormal EBs (panel 2). 
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Fig. 19. — The mass-mass relation for the detached EBs of the TrES dataset. Each category 
is represented by a different symbol. Note that the ambiguous EBs are plotted twice, where 
only one of the solution can be correct. Notice also that in the equal-component solutions 
are clustered along the diagonal, while the unequal-component solutions are clustered along 
the minimum available mass of the Yonsei-Yale isochrones (OAMq). Some of the ambiguous 
solutions deviate from these clusters due to poor constraints on the secondary eclipse, which 
brings about a large uncertainty. Finally, notice the sparsity of EBs populating the low-mass 
corner of this plot (Mi^2 < 0.75 M©). These systems, whose importance is outlined in §3.H 
were modeled using the Baraffe isochrones. CM Draconis (T-DraO-01363) clearly sets itself 
apart, being the lowest-mass binary in the catalog (circle at bottom-left). 
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Fig. 20. — The relation between the EB discovery yield (the fraction of LCs found to be 
EBs) and the absolute value of the Galactic latitude, or for the 10 TrES fields used in 
this catalog (see Tables [1] and [2]) . The solid line is the linear regression of the log of the EB 
discovery yield (r^ = 0.867). Some of the residual scatter can be explained as being due 
to differences in the duration of observations in each field. By including the duration in a 
bi-linear regression, we get a substantially improved fit (r^ = 0.911). 
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Fig. 21. — Low mass candidates (Mi < O.75M0), with their best-fit MECI models (sohd 
hne) . 
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Fig. 22. — The period-eccentricity relation. The lower end of the error bars were truncated, 
where needed, by the measured lower limit, |ecosci;|. Note the lack of eccentric short period 
systems. The diagonal line is provided to guide the eye. 
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Fig. 23. — The period-primary mass relation for eccentric EBs. We included all systems with 
well-determined masses. The area of the grey circles is proportional to the EB's eccentricity. 
All the curves are theoretical boundaries, assuming that the binary components are both on 
the main sequence and have equal masses (g = 1). The left-most dot-dash line demarcates the 
binary contact limit, and the remaining curves mark systems with increasing circularization 
time (see equation [TO!) . Note the abrupt increase in the circularization time for systems more 
massive than ~ 1.25Mq, at which point the stellar convective envelope becomes radiative, 
and thus far less efficient at tidal dissipation. 



